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Introduction 

This  report  contains  sprae  results  on  finite  deformations 
of  an  ideal  elastic  material. 

The  first  chapter  deals  with  a  purely  geometric  question, 
which  is  quite  independent  of  the  physical  properties  of 
materials.  The  question  is  how  much  a  1-1  transformation  of  a 
domain  D  can  differ  from  a  rigid  motion,  if  the  maximum  relative 
extension  ' 


(1) 


ds  '  -  ds 


ds 


occiorring  is  at  most  e.  The  answer  to  this  question  gives 
information  on  the  size  of  the  strains  that  are  bound  to  occ\ir 
in  a  solid  \-jhen   the  solid  is  subjected  to  a  deformation  that 
changes  the  distance  of  one  pair  of  points  by  a  certain  amount. 
It  is  obvious  that  for  a  solid  of  diameter  /\^   a  deformation 
with  maximum  extension  e  cannot  increase  the  distance  of  any 
two  points  by  more  than  e/\  ,  It  is  not  so  obvious  that  mutual 
distances  between  two  points  could  not  decrease  by  larger 
amounts,  IJliat  distance  changes  are  compatible  with  extensions 
of  maximum  size  e  depends  on  the  shape  of  the  body.  The 
geometric  analysis  is  carried  through  here  only  for  the  case  of 
a  square  plate  of  side  a  and  thickness  h,  where  h  <  ^«  ^^   is 
proved  t)iat  there  exists  a  constant  M  with  the  following 
property.  For  any  transformation  x'  =  f(x)  of  the  plate  (in 
vector  notation)  there  exists  a  rigid  motion  x*  =  g(x)  [g(x)  is 
a  linear  transformation  with  constant  coefficients]  such  that 

(2)  |f(x)  -g(x)|  <  m&^/h     . 

The  best  value  for  M  is  not  determined  here,  ,  It  depends 
doubtlessly  on  the  range  of  e-values  admitted.   If,  say,  only 
values  of  e  below  l/lOO  are  considered,  better  estimates  for  M 
6an  b,e  obtained,  than  if  we  permit  for  e  all  values  between  0 
and  1,  For  practical  purposes  it  would  seem  worthwhile  to 
obtain  the  best  M  for  the  limiting  case  e  — >0, 
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Formula  (2)  shows  that  a  deformation  of  a  plate  differs 
from  a  rigid  motion  by  a  quantity  of  order  ea,  unless  the 
thickness  h  of  the  plate  is  small  in  comparison  to  its  width  a. 
It  is  shown  for  thin  plates  the  components  of  displacement 

tangential  to  the  plate  stay  of  order  ea  (at  least  for 

2  2  2 

e  «  h  /a  ) ,  while  the  normal  components  can  be  of  order  ea  /h. 

The  proof  of  (2)  is  based  on  the  realization  that  there 
exist  subsets  called  "cores"  of  the  domain  D  which  is  under- 
going the  transformation,  in  which  the  j?elative  change  in 
distance  of  any  tx^o  points  is  at  most  e.  In  fact,  any  convex 
subset  of  D  whose  distance  from  the  boun.dary  of  D  exceeds  a 
certain  amount  6  =  5(e,D)  is  such  a  core.  Moreover,  5  tends 
to  0  for  e  — *►  0,  This  implies  that  for  convex  D  and  small  e 
there  are  core^  of  D  which  fill  all  of  D  except  for  a  thin 
boundary  layer.   Since  in  a  core  the  mutual  distance  between 
any  two  points  changes  by  a  small  amount  for  small  e,  the  trans- 
formation of  a  core  is  essentially  rigid,  which  leads  to  the 
desired  result.  Similar  estimates  can  doubtlessly  b,e  obtained 
for  D  which  have  shapes  differing  from  square  plates. 

We  have  here  measured  the  amount  by  which  the  trans- 
formation f  differs  from  a  rigid  motion  g  by  the  maximum  of  f-g. 
It  is  more  customary  to  measure  the  difference  by  the  difference 
of  the  matrices  of  the  first  derivatives  of  f  and  g.  Let  p 
denote  the  matrix  formed  from  all  first  derivatives  of  the 
components  of  f  ("Jacobian  matrix").  The  matrix  p  can  be 
decomposed  into  an  orthogonal  matrix  c  giving  the  local  rotation 
and  a  symmetric  matrix  l+v,  determining  the  strain: 

cp  =  l  +  ki   , 

For  small  elongations  y^  is  small.  For  vq  =  o  we  have  a  rigid 
motion  and  c  must  be  a  constant  orthogonal  matrix.   The  question, 
then  is  how  much  c  can  differ  from  a  constant  matrix  for  small  vo, 
N     It  is  shown  by  a  counter  example  that  small  extensions  do 
not  guarantee  that  the  local  rota.tions  at  all  points  are 
approximately  the  same  everywhere.  However,  this  will  be  true 
in  the  mean.  We  can  show  that  in  the  case  of  a  square  plate 
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subjected  to  a  transformation  with  maxim-um  extension  e  there 

is  a  constant  orthogonal  matrix  y  such  that  we  have  the  estimate 


(c  -y)(c  -r)'"'dV  =  0(sa^/h) 


for  the  norm  of  c-y.  This  means  that  on  the  average,  p  differs 

little  from  a  constant  orthogonal  matrix  for  small  e.  For  h  of 

3  1 
the  same  order  as  a  the  norm  is  of  order  ea-'. 

Chapter  II  gives  an  exposition  of  the  classical  equations 

for  ideal  elastic  isotropic  materjLals.  The  relations  are 

equivalent  to  tho^e  given  by  P,  D,  Murnaghan  [2],  Novozliilov  [1] 

and  by  Rivlin  [6],  The  emphasis  is  however  on  arriving  at  a 

formulation  that  permits  formulation  anA  solution  of  problems 

piirely  in  terms  of  Lagrange  coordinates.  The  strains  and 

stresses  which  are  the  quantities  usually  in  the  foreground 

are  replaced  here  by  their  more  immediate  Lagranglan  analogues 

p   =?^      a   =9liM 
^ik   ax^^  '         ^ik   5p.j^ 

where  U  is  the  strain  energy  function.  The  advantage  of  working 
with  these  quantities  (not  sufficiently  stressed  in  the  litera- 
txire)  is  that  the  equilibrium  equations  and  compatibility 
conditions  are  the  simple  linear  relations 

^ik,k  =  °  »    Pik,r  "^  Pir,k  ' 

and  that  the  boundary  conditions  are  all  takpn  on  the  known 
boundary  of  the  body  in  the  unstrained  state.   No  curvature 
tensors  have  to  be  considered,  and  tensor  calculvis  can  be 
avoided  entirely  in  favor  of  matrix  computations.  Non-linearity 
enters  pnly  through  the  formulae  expressing  the  q.,  in  terms  of 
the  p 


rs' 
The  relations 
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between  p  and  q  which  take  the  place  of  the  stress  s.train 
relations  are  analyzed  in  more  detail  in  Chapter  III,   In 

For  a  discussion  of  orders  of  magnitudes  of  rotations  and 
extensions  see  Novozhilov  [1], 


particular,  the  expansion  of  q  in  terms  of  p  near  p  =  1  is 
carried  out  up  to  terms  of  third  order  in  p-1  since  information 
on  these  terms  is  needed  for  the  later  discussion  of  buckling 

of  plates,- 

In  Chapter  IV  we  consider  changes  in  potential  energy 
resulting  from  changes  in  the  transformation  x'  =  f (x),  ,  The 
potential  energy  is  stationary  for  the  equilibrium  state.   It 
is  sho^wn  that  the  equilibrium  state  actually  corresponds  to  a 
minimum  of  potential  energy  if  displacements  are  prescribed  on 
the  bovmdary,  provided  that  for  the  deformations  admitted  the 
matrix  p  is,  restricted  to  a  small  finite  neighborhood  of  the 
unit  matrix,  ThJ.s  yields  a  uniqueness  theorem  for  the  problem 
of  determining  the  strained  state  of  the  body  for  given  dis- 
placements on  the  boTondary  for  an  ai-'bitrary  choice  of  strain 
energy  function. 

Chapters  V  and  VI  derive  th,e  equations  for  bending  of  thin 
plates  using  Lagrange  coordinates.  In  essence  the  method  is 
similar  to  that  used  by  Chien  and  Synge  (see  [3],  Ik-]    )• 
The  plate  in  the  vinstrained  state  is  taken  to.  be  of  constant 
thickness  2h  x^ith  its  faces  given  by  x^  =  i  h.   It  is  assumed 
that  the  transformation  corresponding  to  the  strained  state 
can  be  expanded  into  a  formal  power  series  in  terms  of  x_  and  h. 
We  obtain  from  the  equilibriimi  equations,  compatibility  condi- 
tions and  boundary  conditions  on  the  faces  recursion  formulae 
for  the  cpefficients  of  the  power  series  which  are  functions  of 
x^  and  Xo,  These  reciorsion  formulae  essentially  determine 
successively  the  coefficients  and  excludp  p.   buckled  solution 
unless  on  the  limiting  middle  siirface,  i,e,  for  x-.  =  h  =  0  the 
determinant 

a^u(p) 


3Pik  ^Prs 

vanishes.   The  determinant  vanishes  in  particular  when  the 
stresses  vanish  to  lowest  order  in  h  on  the  middle  sxirface 
x_  =  0,  In  that  case  the  solution  of  the  recvirsion  formulae 
follows  the  pattern  familiar  from  other  asyTiptotic  expansion 
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problems.  At  each  stage  of  the  expansion  the  recTirsion  formulae 
determine  the  expansion  coefficients  only  within  an  arbitrary 
fxmction;  differential  equations  for  that  arbitrary  function 
are  only  found  by  considering  the  reciu?sion  formulae  for  the 
next  higher  orders  and  expressing  that  those  are  compatible 
with  existence  of  a  solution.  In  particular,  the  complete  ' 
system  of  differential  equations  for  the  first  non-vanishing 
quantities  constitutes  then  what  may  be  considered  as  the 
differential  equations  for  thin  plates. 

In  Chapter  V  this  program  is  carried  out  for  the  case  of 
two  dimensional  transformations,  admitting  large  local  rotations 
for  h  — >  0,  The  resulting  differential  equation  for  the  inclina- 
tion 9  of  the  limiting  middle  surface  is  similar  to  the  equation 
of  the  elastica.  In  Chapter  VI  three  dimensional  transforma- 
tions are  considered,  in  which  however  the  middle  surface  shall 
in  the  limit  for  h  — >  0  coincide  with  the  original  unstrained 
surface  x-  =  0  (instpad  of  coinciding  with  a  more  general 
developable  surface).  The  fiirther  assumption  that  the 
t^gential  stresses  are  of  order  h  leads  to  the  equations  of 
v,Karman  and  P8ppl,   It  turns  out  that  the  differential  equa- 
tions for  the  lowest  order  non-vanishing  quantities  only  involve 
two  material  constants,  the  Lame  constants  X  and  [i.     In  the 
derivation  of  the  equations,  higher  order  tjsrms  in  the  expansion 
of  q  in  powers  of  p-1  have  to  be  considered.  These  terms 
involve  in  addition  to  \   and  [i   seven  further  constants,  which, 
however  for  reasons  of  symmetry  do  not  enter  the  final  result. 


Chapter  I 
Geometric  Analysis  of  Deformations  with  Small  Strain 
Let  a  transformation  be  described  by  the  equations 

(1.1)  y^l   =  f^(x^,X2,x^)  ,         1  =  1,2,3, 
or  in  vector  notation  by 

(1.2)  x'  =  f (x)   . 

We  associate  with  the  transformation  the  matrix 

<^-3>  P=(Pik'=(a4) 

and  assume  that  p  has  a  positive  determinant  |p|.  The  corres- 
ponding local  linear  approximation  to  the  transformation  is 
then  given  by 

{l,k)  dx*  =  pdx  , 

"^^^  elongation  of  a  line  element  is  represented  by  the 
expression 


^s    si       dx'-"dx 

where  the  asterisk  denotes  transposition. 

The  elongation  in  different  directions  is  completely 
determined  by  the  positive  symmetric  matrix 

(1.6)  g  =  p'"p  • 

For  dx  varying  in  direction  e  is,  stationary  when  dx  has  the 
direction  of  an  eigenvector  of  g»  The  stationary  values  of  e 
are  given  by 

(1.7)  yT-i 

where  X  is  an  eigenvalue  of  g.  Following  K,  0,  Priedrichs  we 
can  introduce  as  strain  matrix  the  matrix  y|  determined  by  the 
conditions  that  l  +  v^    is  positive  syiranetrlc  and  that 
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(1.8)  ,  (1  ^\)^   =  g  =  p-=^p  . 

The  principal  elongations  e  ai^e  just  the  eigenvalues  of  w  •  If 
we  introduce  the  matrix  c  by 

(1.9)  cp  =  1  +  1^ 

we  have 

p''p  =  (l+>^)^  =  (l+^|■"')(l  +  ^?)  =  p'^c-"cp 
and  hence 

c^'c  =  1  , 

The  matrix  c  ^^d.ll  be  orthogonal  ajid  of  determin,ant  +1,  since  p 
and  1+h  have  positive  determinant.  Equation  (1,9)  corresponds 
to  the  decomposition  of  the  infinitesimal  transformation  into  a 
rotation  c"''  and  a  pure  dilatation  represented  by  1+fe: 

p  =  c''-(l  +yj)   . 

If  the  transformation  is  isometric  we  have 

g  =  1  ,   e  =  0  ,    V|  =  0  ,   p  =  c"'  , 

In  that  case  f  must  represent  a  rigid  motion  and  the  local 
rotation  must  be  given  by  a  matrix  c"  with  constajit  elements. 
The  local  rotation  is  then  the  same  at  all  points.  In  this 
chapter  we  shall  investigate  what  happens  if  one  assumes  that 
the  elongation  e  is  small  for  all  line  elements.  We  shall 
estimate  how  m.uch  the  transformation  f  can  differ  from  a  rigid 
motion  in  this  case, 

Asst:irae  then  that  we  apply  a  transformation  x'  =  f  (x)  id.th 
f  of  class  C  to  the  points  of  a  region  R  in  x-space.  Let  there 
be  given  a  fixed  quantity  e  such  that  for  all  x  in  R 

(1.10)  -E  <  e  =  ^^LzM.  <^   +B 

where  0  <  e  <  1, 

To  illustrate  the  behavior  to  be  expected  we  consider  two 
examples  of  transformations.  The  first  example  to  be  discussed 
will  show  that  even  though  the  elongation  is  small  everywhere 
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the  local  rotation  does  not  have  to  be  approximately  the  same 
In  all  points.  We  nake  use  of  cylindrical  coordinates  r,  ©,  zj 

X,  =  r  cos  Q  ,   X2  =  r  sin  9  ,   x^  =  z 

and  consider  a  transformation  of  the  form 

(1.11)       r'  =  r  ,   ©*  =i  e-Piv)    ,   z'  =  z 

with  a  suitable  function  F(r),   The  elements  of  length  are  giVeh 

by 

2        222      2 

ds     =  dr    +r  d©    +dz'^ 


Here 


ds'^  =  ds^+r^(de'    -  d©^)      . 

|r^(de'^ -d©^)  I   =    jrF'(-2rd©dr  +rP'dr^)| 

<   |rFM(l+  IrP'hds^  <  eds^ 


if 


(1.12)  IrP'(r)|    <  I       . 

It  follows   that  the  elongation  e  satisfies    (1,10),   if   (1,12) 
holds. 

The  matrix  p  is  here  given  by 

/  cos  P+rP'  sin  ©'  cos  ©    sin  P+rP'  sin  ©'  sin  ©    0 

p  =  I  -sin  P-rP*  cos  ©'  cos  ©    cos  P-rP'  cos  ©'  sin  ©    0 

V        0  0  1 

Taking  for  c  the  orthogonal  matrix 

/^-  cos  0  -sin  0         0 

c  =   +sin  0  cos  0         0 

V   0         0      1. 

the  matrix  cp  will  be  symmetric,  when 


I 


^"rl 


^  =  F(r)  +arc  tan  (|  rP'(r))  . 

If  we  now  choose  for  P(r)  the  fimctlon 

F(r)  =  e  log  (log  i)   , 

condition  (1,12)  will  be  satisfied  for  r  <  e"  ,  while 

lim  ^  =  (3D  , 
r— >  0 

so  that  the  rotation  angle  ^  takes  every  value  in  a  neighbor- 
hood of  the  origin. 

This  shows  that  arbitrarily  large  variations  in  local 
rotation  are  compatible  with  small  strains. 

The  next  example  indicates  the  magnitudes  of  the  dis- 
placements in  various  directions  to  be  expected  in  a  trans- 
formation of  ,a  rectangular  plate  when  the  maximum  elongation 
is  of  order  e.   Let  the  plate  be  given  by 

(1,13)    0<x,  <a,   0  <  ^2  -  ^   *        ^~^-il^~^  • 
We  apply  the  transformation 

xJ  =  (R-x^)sin(x,/R)  ,   X2  =  X2  ,   x *  =  R  -  (R-x-)cos (x,/R) 

which  takes  the  plate  into  a  portion  of  a  cylindrical  shell. 
Here  the  parameter  R  shall  have  a  value  exceeding  h  and  a.  The 
strain  matrix  has  the  simple  form 

0    0  \ 
vo  =  I    0      1    0 

0      ij 

in  this  case,  so  that  the  maximum  elongation  is  e  =  h/R,  The 
maximum  horizontal  displacement  occurs  for  x,  =  a,  x_  =  h  and 
has  the  value 

(1,12a)   max  |x|-x^|  =  a  -  (R-h)sin(a/R)  =  ae(l+0(a^e/}i^) )   . 

The  maximum  vertical  displacement  occurs  for  x^  =  a,  x.  =  0  and 
has  the  value 


10 

(1.12b)    max  |x*-x-|  =  2R  sln^(a/2R)  =  0(a^e/h)   . 

Since  here  by  assumption  ae/h  =  a/R  <  1  the  maxirauin  total  dis- 
placement  of  any  point  is  of  order  a  e/h.  In  particular,  if  h 

is  of  the  sarae  order  as  a  the  maximum  displacement  is  of 

2   2 
order  ea.  If  h  is  small  compared  to  a  but  a  e/h  is  still 

bounded  the  horizontal  displacement  is  still  of  order  ea, 

whereas  the  vertical  displacement  goes  up  to  the  order  (a/h)ea. 

This  latter  situation  arises  when  the  total  displacemefit  of 

any  point  does  not  exceed  the  thickness  h  of  the  plate. 

We  find  then  in  this  example  that  for  a  rectangular  plate 
with  sides  of  the  same  order  a  and  maximum  elongation  e  the 
displacement  is  of  order  ea.  The  horizontal  displacement  is 
still  of  the  same  order  as  long  as  ^e  does  not  exceed  the  ratio 
of  thickness  h  to  width  a.   In  any  case  the  total  displacement 
is  of  order  ea  /h,  as  long  as  e  does  not  exceed  h/a.  We  shall 
find  essentially  the  same  result  for  general  transformations  of 
the  plate  with  raaximvim  elongation  e  (af,ter  applying  a  suitable 
rigid  motion  to  the  transformed  region) » 

Let  then 

X*  =  f (x) 

be  a  transformation  of  a  region  R  with  boundary  B,  where 

for  x€R,  If  C  is  any  curve  in  R  and  C'  its  image  under  f,  we 
have  the  inequalities 


(1  -  e)  r  ds  <  r  ds '  <  (1  +  e)  j 
C      C»  C 


ds 


for  the  lengths  of  the  curves  C  and  C', 

Let  P,  Q,  be  any  two  points  of  R  such  that  the  line  segment 
C  with  endpoints  P,  Q,  belongs  to  R,   Let  PQ  denote  the  distance 
of  P  and  Q,   Then 

P'Q'  <  J  ds'  <  (1+  e)  j  ds  =  (1+  e)  PQ  . 
C»  C 


K'Kf    ^]-. 
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If  in  addition  the  line  segment  y'  with  endpoints  P',  Q*  belongs 
to  the  image  region  R',  then  y'  is  the  image  of  an  arc  y 
joining  P  and  Q,  and 


r 


P»q'  =    ds'  >  (1-  e)   ds  >  (1  -  e)  PQ 
Hence 


Y'  Y 


(li.3)         (1-  e)PQ  <  P'Q'  <  (l  +  e)PQ  , 

if  P  and  Q  are  such  that  the  segment  PQ  belongs  to  R,  the 
segment  P'Q.'  to  r'. 

Consider  now  the  case  of  a  line  segment  PQ  lying  in  R 
which  is  such  that  the  segment  P'Q'  does  not  lie  in  R,*  Let  S 
be  a  point  on  PQ,  such  that  P's'  belongs  to  R'.   The  set  of 
these  points  is  closed;  it  also  contains  P  but  not  Q,   There 
is  then  a  point  S  in  this  set  furthest  away  from  P,   The  seg- 
ment p's'  belongs  then  to  r'  but  cannot  lie  completely  in  the 
interior  of  R',   Thus  P'sJ  contains  a  point  T'  of  the  boundary 
B'  of  R'.   T'  is  the  image  of  a  point  T  on  the  boundary  B  of  R, 
Let  y  denote  the  original  in  R  of  the  segment  y'  joining  P* 
and  S',   Since  y  passes  through  T  we  have  (see  Figure  1) 


(l-e)(PT +T3^)  <  (1-e)  |  ds  <  f  ds »  =  p's^  <  (l+e)PS^   . 

Y      Y' 

Then  also  ^ 

(l-e)(PT+ TQ)  <  (l-s)(PT+ TS^  +  S^)  <  (l+e)PS^  +  (l-e)S~Q 
<  (l+e)(PS^ +S^)  =  (l+e)PQ  , 
Hence  there  must  exist  a  point  T  on  the  boundary  B  such  that 

PT +TQ  <  ^  PQ  . 

The  points  T  satisfying  this  inequality  fill  a  certain  ellipse 
with  foci  P,  Q  and  minor  semi-axis 


a.l.i- 
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This  Is  also  the  largest  distance  any  point  P  of  the 
ellipse  can  have  from  the  segment  PQ,  It  follows  that  for  a 
segment  PQ  In  R  the  Inequality  (11)  holds,  If  every  point  of 
the  segment  PQ  has  a  distance  from  the  boundary  B  which  exceeds 


1-e 


PQ 


Pig.  1 


We  call  a  subregion  D  of  R  a  core,  if  D  is  convex  and  If 
the  convex  hull  of  the  Image  D'  of  D  Is  contained  in  R*.   If  D 
is  a  core  then  the  inequalities  (11)  hold  for  any  two  points 
of  D,  We  have  from  the  preceding  argument  that  any  convex  sub- 
set D  of  R  is  a  core  if  the  minimum  distance  of  D  from  the 
bovindary  B  of  R  is  at  least 

1-e  •^  » 

where  /\  is  the  diameter  of  D, 

It  will  be  seen  that  the  transformation  f,  can  be  approxi- 
mated by  a  rigid  motion  in  the  points  of  a  core.  Let  there  be 
given  a  core  that  contains  k   points  A  ,  A, ,  Ap,  A-  where,  the 
directionja  A  A^  for  1  =  1,2,3  are  mutually  perpendicular.  Let 

a.  =  A~a7, 
1    0  1* 


We  make  the  assumption  that 


-.< 


-  -vV,'.   7  >^ 
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(l,li|.)  a,  =  a2  >  a,  >  0  , 

In  a  suitable  rectangular  coordinate  system  the  A.  will  have 
coordinates 

A^=  (0,0,0),   A^=  {a^,0,0),   A2=(0,a2,0),   A2=(0,0,a^)  , 

By  a  suitable  orthogonal  transformation  applied  to  R'  we  can 
bring  about  that 

A^=  (0,0,0),  A^=  (a^^,0,0),  A^  =  (a2^,a22,0) ,  k^^  =  {a,^^,&^^,ei^^) 

where  a..  ,   Sipp'   ^-^-^  ^^®  positive. 

For  any  point  P  of  the  core  we  put 


^i  ~  ''^i^  '   ^i  "^  ^i^'        ^^^  ^  ~  0'^»2»3  • 
Then  for  P  =  {x^,x^,x^) ,   P'  =  (x^,x^,xj)  ,  i  =  1,2,3 

2^i^i  =  ^o  -  ^i  "■  ^i  '   ^Vl   =  ^f  -  ^i^  "■  4 

and  hence 

(r'  -r  )(r'  +  r  )  -  (r'-r.  )(r.'+r,)  ■ 
(1.15)  4'^,=      "      "      "    2I1 —  '   '   '   . 

Here  for  i  =  0,1,2,3 


(1.16)   r^  =  ApT  =  A^P'  +  0(A^A^)  =  A^P  +  0(eA^P)  +0(A^A^) 


=  r^  +0(er^)  +0(A^A^)   , 
In  particular  for  i  =  0 

(1.17)  r^  =  r^+0(er^)   . 

If  we  take  here  for  P  the  point  A-  relation  (1,16) 


becomes 
and  hence 


a^^  =  a,  +  0(sa,  ) 


Vl  =  '^11-  ^^l'  =  o^'^^i^  • 


fl 


l     '?.o  ^ 


1^ 

We  restrict  P  to  points  of  the  core  with  r  <  i|.a, ,  Then  also 


r.  =  A.  P  <  A.A  +A^P  =  a.  +  r^  <  ^a^      , 
1     1   —   lO    o      1    o—    1 

We  have  then  from  (1,16)  for  1=1 

r^  =  r^  +  0 ( er,  )  +  0 ( ea^ )  =  r^  +  0 ( ea^ )  =  0 ( a^ ) 

It  follows  from  (1,15)  that 

(1.18)  x|  -Xj^  =  0(ea3^)   . 
In  particular  for  P  =  A2,A 

(1.19)  &2i   ~  O(ea^)  ,    a   =  O(ea^)   , 

Taking  P  =  A2  we  have  from  (1.17),  (1.19) 


^22  ^v/^22  "^  ^21  "^  ^^^21^  ~  rQ  +  0(ea^) 
=  r  +  0  ( ea^ )  =  02  +  0  ( ea^ ) 


and  hence 


V2=/4l-'^^22-^2^^  =  0(ea3^)   . 


Then  generally  for  points  P  of  the  core  with  r  <  i;a,  by  (1,16) 
and  (1,15)  (using  a^   =  o.-^) 


^2  ~  ^2+0(80^)  =  O(a^) 
(1.20)         X2  -  X2  =  O(ea^) 

In  particular  for  P  =  A-  relation  (1.20)  becomes 

a^2  =  O(ea^)   , 
and  hence  also 

^33  "- J "-33  "  ^31  "■  ^%  "   °(^^31  -^  ^32  J  =  K"  °^^^1> 
=  ro+ O(ea^)  =  a^  +  o(ea^)   , 


•■\     AC!,:     ■' 


•TOO    ar. 


'v..f 


'■  ..  ■  .)■ 


r-  1 


.  fX  ,  •' ) 
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It  follows  that 


A^A^  =y  a^]^ +a^2  •*'^^33 -^3^^  =  O^^'^l^   • 
Then  from  (1,15),  (1.16)  for  any  P  of  the  core  with  r^  <  l+a^^ 

rl   =  r^  +0(ea^)  =  O(a^) 

(1.21)  X-  -  x^  =  0(ea,/a-)    , 

Take  now  the  case  where  R  Is  the  rectangular  plate 
(1,21a)        0  <  x^  <  a.   ,       0<Xp<a,   0<x^<h 

with  0  <  h  <  a.  For  given  positive  5  the  cell  R^ 

(1.22)  5  <  X,  <  a-5  ,    5  <  X2  <  a-6  ,   6  <  ^^  <  h-6 

will  have  distance  6,  from  the  boundary  B  of  R,  and  its  diameter 
villi   be  at  most, /3a,   It  will  be  a  core,  if 

Ass^jme  that  e  is  so  small  that 

(1.23)  ^^<p(l-^)^ 
and  take 

(1.2i|)  S  =^  a  . 

Then  5  <  h/I;.,   The  cell  Rg  has  sides 

a,  =  a  -  25  ,  ap  =  a  -  25  ,   a-  =  h  -  25  , 

and  hence  is  not  empty.  Moreover 

cuj_  =  0(a)   ,  a^/a^  =  0(a^/h)   , 

It  follows  from  (1,18),  (1,20),  (1,21)  that  after  a  suitable 
orthogonal  transformation  applied  to  R*  we  have  at  least  for 
the  points  of  the  core  R. 


16 

(1,25)   x^-x^  =  0(ea)  ,  x^-Xg  =  0(ea)  ,  x^  -  x^  =  0(ea  /li)  , 

These  are  exactly  the  orders  of  magnitude  for  horizontal  and 
vertical  displacements  that  occurred  in  the  special  example 
(1.12  a,b)  when  ea  /h  was  bounded.   They  are  proved  here  only 
for  the  core,  whic)i  however  in  the  limit  for  e  — >  0  will  fill 
the  whole  region  R, 

A  different  kind  of  argument  is  needed  to  estimate  the 
displacements  outside  the  core  R_  and  also  to  tj?eat  the  case 
where  h  is  so  small  that  (1.23)  no  longer  holds.  For  this 
piArpose  vie   consider  first  the  case  of  a  core  which  contains  a 
cube  of  side  a  and  vertices  A  ,A_,,,,  ,,  If  P  is  a  point  of  the 
core  with  A^  =  a  we  have  by  (1,18),  (1,20),  (1.21)  after  a 
suitable  rigid  motion  applied  to  R* 

PP^  =  0(ea)  , 

If  Q   denotes  the  angle  PA  P'  we  have 

^   PA^P'  =  fi  =0(PPVAoP)  =  0(e)  , 

If  Q  is  the  point  of  distance  a  from  A  on  the  ray  from  A 
opposite  to  P  and  if  Q,  also  belongs  to  the  core,  then 


7-  'o(^_'oC 

and  hence 


^  QAqQ»  =  0(e) 


^  P'A^Q'  =  ^   P»A^Q»  =  71+ 0(e)   . 

Thus  if  PQ  is  a  line  segment  of  length  2  a  and  A  its  midpoint 
we  have 

^  P»A^A»  =  7t+0(e)   , 

provided  there  exists  a  core  that  contains  P,  Q  and  also 
contains  a  cube  of  side  a  and  vertex  A  ,  This  is  true  regard- 
less of  the  special  rigid  motion  applied  to  R'  which  would  not 
change  the  angle^,.  , 

Let  Po*^l*****^n  ^^  equidistant  points  on  a  line  segment 

L  =  P  P  of  R,  where 
on* 


^i^i+l  =  ^   • 


^■v-  ^.i?-. 


,•    '•>. 


.     i2)J 


1   tr 


uJ 


17 

Let  there  exist  for  each  1  =  l,2,,,,,n-l  a  core  containing 

^i<-l  ^"^  ^i+1  ^^^  also  containing  a  cube  of  side  a  and  vertex 
P.,  Let  simi;Larly  P  ,  P,  belong  to  a  cubical  core  of  side  a 
and  vertex  P  ,  The  angle  between  two  successive  vectors 


PjPi^n  is  then  0(e),  Hence  the  angle  between  any  vector  P.'pJ  , 

and  P'P^  is  0(ne),  We  assumed  that  there  was  a  cubical  core 

of  side  a  and  verte;x  P  ,  Let  A,  ,  A^,  A_  be  the  vertices  of  the 

core  adjacent  to  P  ,  We  take  P  as  origin  and  the  adjacent 

edges  of  the  cube  as  coordinate  axes.   By  a  suitable  rigid 

motion  we  bring  about  that  P*  =  P  ,  A*  lies  on  P  A^ ,  a'  lies  in 

,  o    o'   1         o  1'   2 

the  plane  P  A.A,,   Then  the  angle  between  P  P»  and  T^T  ,  that 

O    i.    c.  .ox  OJ.' 

-> 

is  the  angle  betv/een  P'P»  and  L  is  0(e),  It  follows  that  the 

>   o  1 >     -.— > 

angle  between  P'P«   and  L  or  between  P'P'  ,  and  P.P. .,  is 
i  i+1  i  1+1      i  1+1 

0(ne),  Since  also 


^Fi+i  =  Vi+1 -^  0^^^^     ' 
we  have  for  P^  =   (x^^,x^2»^i3) »   ^i  ~   ^^il'^12*^i3^ 

^1+1  k  -  ^ik  =  ^1+lk-  ^ik-'^^"'^^     ^°^     ^  =  l'^'^   • 


ik 
Hence 


^i+1  k  -  -Ik^   =  2^   ^^i+1  k  -  ^ik^  +0(n2ea) 


=  x^j^  +0(n^ea)   , 
Consequently 

Fp^  =  ©(n^ea)   . 
Let  again  R  be  the  region 

0<x^<a,        0<X2<a,        0<x^<h<a     , 
Let  Rq     be  the  subset  of  R  described  by 

^1  <  ^1  <  a-^i    »        5]^  <  X2  <  a-5^    ,        6^  <  X-   <  h-5^      , 
where 


•  one. 


Then 
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t 

5^=  ye  h    , 

5,    <  h/lj   if  we  assiHTie   that   e  <  l/l6.  Take  for  P     the  point 


(5,, 5^, 5^)  and  for  Q  any  point  of  Rg  •  We  divide  the  line 

segment  L  =,  P  Q,  into  n  equal  parts  a  by  points  ^ \>^ 2* "  *  *^ n-l 
with  P  =  Q,   V/e  take  here  for  n  some  integer  satisfying 


20  g  <  n  <  20  1  +  1  , 


Then 


^  =  S^o^^l^^^^^^l^^-2\^   • 


Every  point  P.  of  L  is  the  vertex  of  a  cube  in  Re   v;i,th  sides 

of  length  a  which  are  parallel  to  the  coordinate  axes,  ,  The 

convex  hull  of  this  cube  and  of  P.  ^    and  P. .t  is  a  core.  For 

i-1      i+1 

its  diameter  /\^   is  at  most 

(l+yi)a  <^   (l+yi)h 
and  its  distance  from  B  at  least  5. ,  so  that 

Consequently  (after  the  proper  orthogonal  transformation  applied 
to  R') 

QQ'  =  O(n^ea)  =  O(nep^)  =  O(nea)  =  0(ea^/h) 

for  any  point  Q  of  R^  . 

°1 
It  remains  to  extend  the  inequality 


QQ*  =  0(ea2/h) 


to  all  points  of  R, 
that 


If  Q  is  any  point  of  R  there  exists  a  point  P  in  R.  such 


PCi  <  5   , 
Then 


Q*Q  <  Q'P*  +  P'P  +  PQ  <  (1+  e)PQ  +P*P  +PQ  <  (2  +  e)5  +P'P  , 


»    ' '  ^.  "  J^'- 


•I  -.  • !  ■ 


U    \ 


■■} 


)    T-C^"--' 
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If  here  the  inequality 

p'p  =  0(ea^/h) 

has  been  established  already  for  the  points  P  of  R„,  and  if 
6  =  0(ea^/h)  it  follows  that 

qTq,   =  0(ea^/h) 

for  the  arbitrary  point  Q,  of  R, 
Taking  here 

5  =  5^  =  /e  h 

we  find  that  (1,26)  holds  generally  in  R,  if 

ye  h  =  0(ea  /h)     , 

This  is  the  case  when 

/e  — w  >  1  , 
h  ■" 

Take  on  the  other  hand  the  case  where 

„2 

(1.27)  •  v^~-5  <  1  . 

h"^  "■ 

For  every  Q  in  R  there  exists  a  cube  C  of  side  26^  wh^-ch  has  Q 
as  a  vertex  and  sides  parallel  to  the  coordinate  axes.  Let  the 
distance  6^  of  Q,  from  the  boundary  B  of  R  satisfy 

The  vertex  A  of  C  opposite  to  Q  will  have  distance  at  least' 

25^  from  B  and  will  lie  Inside  R.  ,   In  fact  the  subcube  C.  of 
1  o,  ,         1 

side  6-  and  vertex  A  will  lie  completely  in  R.  , 

dm  O  O^ 


J      ,eoXi}   0 


,^fO'I  ' 


p 

a^, 

R. 

-^1 

c. 

-.. 

c 

4 

'^/ 

1 

1 

1 

20 


Pig,  2 


Since  C  has  diameter  ^   =  yi2  5^  and  distance  5p  from  B  the 
cube  C  is  a  core  since 


^   A  -  V^^  5   <  5 


Let  A_ ,  Ap,  A  be  the  vertices  of  the  subcube  C.  adjacent  to  A 
Put 


r^  =  A^Q  ,   r^  =  A^Q' 

Then,  since  the  A  belong  to  Rr 
r^  =  0(5^) 


(i  =  0,2,3)  • 


^i  ~  '^i^*  "^  ^1^*  "^  O(A^A^)  =  r^+0(er^)  +0(A^A^) 
=  r^  +  0(e5^)+  0(Ea^/h)  =  r^+0(e^/^h)  +0(ea^/h) 


=  r^+0(ea^^)   , 
It  follows  from  (1,2?)  that 


2   2 


■  ^'  'Q-     ^   =  0(eaV^(l+  ea^'hl^)) 


=  0(eaV^(l+ e^/Vh"^))  =  oCea^A)  . 


Applying  equation  (1.15)  with  a^  =  a2  =  a-  =  5^ 
for  the  coordinates  x.  of  Q  and  xi  of  Q 


21 
we  find  that 


^i  "^i  ~  O(ea^A) 


and  hence  that 


Q'Q  =  0(ea^/h) 


for  any  Q.  in  Re  ,   where 
53 


Since 


^3  -   l^T-  ^1  -   iTT  ^  • 


5-  =  0(£h)  =  0(ea^/h) 


it  follows  then  that  (1,2  6)  holds  for  any  Q,  in  R  (see  page  19). 
Summarizing  we  have  the  following  theorem: 
There  exists  a  constant  M  with  the  following  property: 
For  any  transformation  x'  =  f (x)  with  maximum  elongation  e 
applied  to  a  square  plate  x-jlth  sides  a,  a,  h  (where  h  <  a) 
the  maxiraijim  displacement  of  any  point  after  applying  a  suitable 
orthogonal  transformation  is  at  most 

Mea^/ii  , 

¥e  can  conclude  that  on  the  average  the  matrix  p  will 
differ  little  from  the  \anlt  matrix  and  that  the  lopal  rotation 
matrix  c  will  also  differ  little  from  the  identity. 

First  of  all  the  symmetric  matrix  -q    caji  be  brought  into 
diagonal  form  by  an  orthogonal  transformation.  Since  the 
eigenvalues  of  /•  are  at  most  e  in  absolute  value  and  the 
orthogonal  matrix  has  bounded  elements,  it  follows  that 
r^  =  0(e),  Then 

P  =  c''^(l+  h)  =  0(1+  e)  =  0(1)  , 

4 

Consider  now  any  parallel  to  the  x,-axis.  On  it  ds  =  dx,  and 


I  ds  *   i    2  2     2 

^1  =  Pil^^l  =  Pil<^^  ,    ^r-=JPll  •*•  P21  -^^31 


'    r.      ~. 


-■  ■  ■  * 


IC:  i  '^JO 
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Since  ds V^s  lies  between  1-e  and  1+e  we  have 

(l-e)2  <  p^3_  +  p^^  +  p|^  <  (1+e)^  , 
In  particular, 

I  Pill  <  1+^  • 

Integrating  over  the  portion  of  the  parallel  to  the  x,-axis 
lying  Inside  R  v;e  have 

dx^  =  \(lx^  +  Oi&{a./h))      , 

^  2         ' 

since  the  endpoints  have  moved  at  most  by  amounts  0{e(a  /h)). 

This  relation  can  be  written 

'  2 

I  ^11  ^1  ~  ^"*'  °(^^^  A^^ 

or 

^  4 

I   (l+e  -p^^)dx^  =  0(e(a^^)  . 

Integrating  over  all  parallels  to  the  x,-axis  we  find 

I  /  /  (l+e -Ptt  )dx  =  O(ea^)     (dz  =  dx,  dx^  dx-j)  . 
jJJ  •*"*•  1   <i   J 

R 

Here  the  integrand  is  non-negative.  Since 

|l-p^^|  <  |l+e -P3_3_l  +  |e|  =  (l+e  -p^-_)  +  E 

It  follows  that 

iff   |l-P3_3_|dx  =  0(ea2) 

p.^  being  bounded,  we  have 

(1-P-L3_)^  =  Odl-p^^^l) 
and  thus  also 


II  j    (l-p^]_)^dx  =  O(ea^) 


'R 

We  see  that  in  the  mean  p, ,  differs  little  from  1,  Moreover, 
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P?i  +P?!   5  ^1+^-  Pm)(l+e  +Ptt)   =  0(1+6  -Pn.)     . 


'21     ^31 


ai' 


Hence 


i'l^ 


<pii-i)^-^pii+P3i 


dx  =  O(ea^)      . 


In  the  same  manner  one  finds  that 


ff 


Pl2+(P22-1^^  +  P32 


dx  =  O(ea^)   , 


R   •- 


P^3+P^3-Kp33.1)2 


dx  =  0(e(a^/4i))  , 


Altogether 


I  rr 


!i  I  (p  -l)(p"--l)dx  =  0(e{a^A))   . 


jjj 


We  measvire  the  deviation  of  the  rotation  matrix  c  from 
vinlty  by  the  matrix 

(c-"--l)(c  -1)   . 

Now 

4 

c  =  (l  +  )|)p~^  =  (l  +  i-f^)"^p'"'  =  p-"--  .](1+  y[)'\''  =  p""'+0(e)  . 
Henoe 

j}  (c"''-l)(c-l)dx  =  jll  (p-l)(p'='-l)dx  +  0(ea3)  =  0(e(a^A))  . 


This  shows  that  for  small  e  the  rotatlpn  matrix  c  differs  In 
the  mean  little  from  a  constant  matrix. 


2k. 
Chapter  II 
Equations  of  Motion  and  Equilibrium  In  Lagrange  Coordinates 

A  material,  Is  considered  here  in  two  states,  "strained" 
and  "unstrained".   The  Cartesian  coordinates  of  a  particle  in 
the  unstrained,  state  ("Lagrange  coordinates")  shall  be 
(x,,Xp,Xo)  =  X,  The  Cartesian  coordinates  of  the  same  particle 
in  the  strained  state  ("Euler  coordinates")  at  the  time  t  shall 
be  (x|,x*,x*)  =  x',  where 

x'  =  x»(x,t)   , 

We  introduce  again  the  Jacoblan  matrix 


(2.1)  p  =  (Pii^)        , 

Velocity  and  acceleration  are  respectively  the  vectors  with 
components 


Let  dV  and  p  denote  element  of  volume  and  density,  in  the 
\mstrained  state,  dv'  and  p   that  in  the  strained  state. 
Conservation  of  mass  implies  that  for  the  element  of  mass 

(2.2)  dm  =  pdV  =  p'dV'   , 
wh,ere 

(2.3)  dV»  =  |p|dV  , 

(Here  |p,|  denotes  the  determinant  of  p,  which  is  assumed  to  be 
positive, ) 

The  stress  in  the  strained  state  is  described  by  a  matrix 
T  =  (t^. ,  ),  such  that  the  force  acting  on  a  surface  element  dS* 
with  unit  normal  £;*  in  the  strained  state  is  given  by 

(2.I^)  \\A^^'    • 

In  addition  to  the  stresses  there  may  be  present  external  forces 
represented  by  a  vector  P  =  (F^jEpjE,)*  such  that  PdV*  represents 
the  external  fprce  acting  on  a  volume  element  dV'  in  the 
strained  state. 


-••-■-  / 
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Let  R*  be  a  region  with  jDoimdary  S*  cut  out  from  the 
material  in  the  strained  state^  Let  £,}   be  the  components  of 
the  exterior  unit  normal  of  S*,  Then  Newton's  law  of  motion 
gives 

|i  p»x{dv»  =  J(i  P,dv'  +11  ^,,^^las*    '. 

R*  R*  S» 

Applying  the  divergence  theorem  to  the  surface  integral  and 
sh-rinking  R*  into  a  point  we  find  the  equations  of  motion 

(2.5)  p  x^  =  Pj^  +  -r-T  . 

The  work  done  by  the  surface  stresses  and  exterior  forces 
on  R*  in  unit  time  is 

(2.6)  W=Jj  i^ik^i'^^idS*  +  111  P^x^dv'   . 

s*    '     r'^ 

The  kinetic  energy  of  the  mass  in  R*  is 

(.2.7)  K  =  i!!U;x?p»dv»   . 

R* 

We  now  postulate  the  existence  of  a  strain  energy 
fimction  as  follows: 

Principle:  There  exists  a  function  U  =  U(p)  (the  strain 
energy  per  unit  voliome  of  the  unstrained  state)  such  that 

R 

where,  R  is  the  region  with  boundary  S  in  x-space  corresponding 
to  R*. 

We  have  from  (2,6),  (2,5) 


/      »  .       « 
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(2.9)        W  = 


If 


R 


P 


1        r 


11  ik   ix^   .vxJL 


dT 


JJJ  ''ikP 


PirP^'^dV' 


R» 


where 


p-^  =    (^^^) 


is   the  matrix  reciprocal  to  p.      On  the   other  hand 
M'  I  >.  <  I  I  j 

'ir  ^"^  = 


(2.10) 


_d_ 
dt 


ji  U(P)  dv=  IJ!  ^11 


R 


ir 


""-   '^  =  !J.!  BT  ^ir^ir  ^^'  ' 


R» 


where  the  q.   are  the  functions  of  p  defined  by 


(2.11) 


^ir  -    dP-. 


-vp . 


By  comparison  of  (2.8),  (2.9),  (2.10),  x^e  find  for  R' 
shrinking  into  a  point  uhat 


(2.12) 


.rkj 


1 


'4kP""Pir  -  TpT  ^ir^'lr   ' 


I«Je  raake  the  assumption  that  a  particle,  whatever  its  local 
deformation  p  and  stress  at  a  given  moment,  can  have  any  given 
deformation  rate  n.   Then  (2.12)  implies  that 


rk 


'-ik^- 


TpT 


q^ 


ir 


Solving  these  equations  for  p  we  find 
(2.13)  Xil,  =-l^qirPkr 
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or,  in  matrix  notation, 

(2.i;|.)  t  =  -ri^  qp-::-    . 

Relations  (2. !,'[.)  take  the  place  of  the  stress-strain  relations 
of  the  classical  linear  theory?,  supplementing  the  equations  of 
motion  (2.5).   Vie  next  derive  the  equations  of  motion  and  bound- 
ary conditions  in  Laf,;ranse  coordinates  x. 

VJe  have  frora  (2.10),  (2.3),  (2.6),  (2.7) 

d  ill   .    \\^      •      'l-  .  rff3%x-   • 

^  U  p   dV   =      I    q.    p .      dV  =   \\\  q .    x'.^      dS  -        I  .^-^  xl    dV 

dt  JJ.i     ^  J JJ  ^ir*^ir  jjj  ^ir   i^r  Jii  c^x  i 

R  R  ■  S  R        ^ 

I  (  .  I'i'l'     ^^ip    .     ,     ,-1  dK 

=  Jj   ^ir-l^r  ^S-   |ji    -^xiIpI    IdV    =W-§ 

S  ■  R'        ^ 

O  fij'" 

=  !;  ''^ik^i^i  ^s'  -»■  jlj  [Fi  -  p'xijx!  dv.  . 

S«  R» 

At  any  (^iven  raoment  t  the  x!  are  Independent  of  the  p.,  and  q., 

and  can  be  taken  as  arbitrary  functions.   Taking  first  for  the 

x!  arbitrary  functions  vanishing  on  S  respectively  S'  i-^e 

X  ■  1 

find  by  comparison  of  int6,^rands  that 


(2.15)  p'xl  =  F.  + 


1   ^^ir 


1    !p   'X. 

Comparing  next  the  remaining  surface  integrals  for  xl  arbitrary 
on  S  we  obtain  the  relations 

(2.16)  q.  ^   dS  =  L..  5'  dS'   , 

Here  dS'  is  the  surface  element  in  the  strained  state  with  unit 
normal  £,'  that  originates  from  a  surface  element  dS  with  unit 
normal  E,   in  the  unstrained  state. 

If  U  is  known  in  its  denendence  on  p  and  if  the  P.  are 
knovm  as  functions  of  x  or  x'  and  t,  we  have  in  (2.l5)  a  system 
of  second  order  partial  differential  equations  for  the  x!  as 


These  equations  are  given  by  Novozhilov ^1]  III  (3^),  and 
Reissner  [51  for  the  case  of  time-independent  transforiflations , 
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functions  of  the  x,  .   Stresses  and  strains  in  the  ordinary  sense 

do  not  enter  these  equations  of  motions.   Their  place  is  taken 
by  the  matrices  q  and  p,  which  are  conriected  by  the  relations 

(2  17)  q    =  ^-(P) 

-  ik 

Boundary  conditions  of  the  type  of  restrictions  on  the  T!'.,  on 
S'  yield  by  (2.16)  restrictions  on  the  q.,  on  3. 

In  tlie  special  case  corresponding  to  elastic  equilibrium 
and  in  the  absence  of  external  forces  the  x!  are  independent  of 
t,  and  the  equilibrium  equations  take  the  simple  form 

(2.18)  --^   =  0    . 

oXj^ 

Making  use  of  (2.17)  we  obtain  from  (2.16)  a  quasi-linear  system 
of  first  order  equations  for  the  V^r^i    v/hich  has  to  be  supplement- 
ed by  the  obvious  ''compatibility  conditions"  for  the  p., 

So  . ,    '. n . 

(2.19)  ■  ^'^  -  '  -^^ 


cX       .-..Xi 

r      k 

The  tvjo  most  corainon  ty^Dss  of  boundary  conditions  are  that 
either  the  displacements  or  the  tractions  vanish  on  portions  of 
S.   In  the  first  case  we  have  xi  -x.  ~-   0   or  equivalently 

(2.20)  p,^  -  5.^^  =C7,4k 

with  suitable  proportionality  factors  r. .   In  the  second  case 
"ik^k  ^^'  =  Q  °^  ^y  (2.16) 

(2.21)  q^^£,^   =  0   . 

In  either  case  the  problem  of  determining  the  xl  or  p.,  is 
described  comoletely  in  La.tran'.e  coordinates. 
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Chapter  III 

Invariants  of  the  Strain  Energy  Function 

VJe  assume  that  the  strain  energy  density  U(p)  is  indepen- 
dent of  the  choice  of  the  cartosian  coordinate  system  used. 
This  amoimts  to  the  requirement  that 

(3.1)  u(p)  =  U(cpc-::-) 

for  any  orthogonal  matrix  c.   In  addition  we  restrict  ourselves 
to  an  isotropic  medium.   This  has  as  a  consequence  that  U(p) 
stays  unchanged  if  a  rigid  motion  is  superimposed  on  the  trans- 
for'-aation  x'  =  f(x)  with  Jacobian  matrix  p,  i.e. 

(3.2)  U(p)  =  Ij(cp) 
for  orthogonal  c.   Altogether  then 

\J{pj  )  =  iKp^) 
for  two  matrices  p, ,  Pp  that  are  connected  by  a  relation 

(3.3)  Pg  =  cp^C 

where  c  and  C  are  orthogonal.   The  necessary  and  sufficient 
dondition  for  two  matrices  p..  ,  p^  to  be  connected  by  a  relation 
(3.3)  with  orthogonal  c,C  is  that  the  matrices  P-JP-.  and  p^Pp 
have  the  same  eigenvalues.   Consequently  our  assumptions  imply 
that  U(p)  is  a  synraetric;  function  of  the  eigenvalues  of  p  "p  =  g. 

Any  syraraetric  fuiiction  of  the  eigenvalues  of  a  matrix  g 
can  be  expressed  in  terns  of  any  three  independent  symmetric 
functions  of  the  eigenvalues,   VJe  choose  here  as  basic  symmetric 
functions  tiie 


(3.1;)  s^  =  [(g  -  1)"^]       for  n  =  1,2,3   , 

where  we  denote  generally  by  [a]  the  trace  of  the  matrix  a.   Then 
(3.5)  IJ(p)  =  V(s-j_,S2,S2)   . 
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It  is  easy  to  establish  the  connection  betl^Jeen  the  rela- 
tions given  here  and  tliose  ^iven  customarily,   in  vjhich  the 
strain  matrix  is  t alien  to  be 

K  =  |(s  -  1)       . 

If   X-],?ioAo   denote   the   eigenvalues   of  ^   and   I,,Ip,Io    the   elem- 
entary  syriiinetric   functions   of  those    eigenvalues,    we   have 


/  s. 


(3.5a) 


5  ^\  =  2Ii 


{s^=^k4  =  kill-   2I2) 


\ 


Then  U  =  V(s,  ,Sp,s^)  becomes  a  function  of  the  I,  .   One  can 
write  TJ  in  the  form  U  =  \i/(^)  vjhere  \|/(^)  =  \{f(^'"').   Then 


^HKi  =  'diiiK) 


zKr 


f  or  ^  =  ^ 


'ik         -"^ki 
One   obtains,    using  this   syiranetry  property, 


-^K.  ''^^ 


%-«. 


ou 

Dx 

3P.V 

>     P.w    - 

Ox 

and  from   (2.13) 


^..   =  IpI"'^ 


p-i  .Pi 


'ik        '^■''        ?/^    '  ^is^kr 
^rs 


^-#-  p.    P, 
P    v<^3      is^kr 


It  is  easy  to  find  the  stress  strain  relations  correspond- 
ing to  a  given  choice  of  f'onction  V.   Vie  have  from  (2.17),  (3.5) 


^U  =  '^ik^Pik  =  f^""  ^^^    = 


ds 


n=l 


n 


n 


See   P.    D.    Murnaghan,    [2],    Chapter   III. 


See   Ilurnaghan,    [2],   p.    56. 
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By   (3.1|-) 

ds^  =  nf(G-l)""^   dg]    =  n[(g--l)''-^   d(p"p)] 
=  n[(g-l)^"^    (dp''~")p   +    (s-l)''"^p'''"   (dp)] 
=  2n[(g-l)''"^p"''   (dp)] 

Hence   by   comparison,    for   arbitrary  dp, 

3 
q''  =^  2nVg    (g-l)^^"^p'"* 
n=l  n 

or 

(3.6)  q  =  EI  2nV^  P(p''p-i:"'^   =  XT  2nV^    (pp  •''■-!  )"'^p      . 

n=l  n  n-l  n 

It  follovjs  from  (2.1ij.)  that  the  stress  matrix  t   is  given  by 

3 

(3.7)  t  =  IpI'-'-  ZI  2nV^^  (pp'-'-D'^'V"   . 

n-1    "n 

This  formula  makes  evident  that  the  stress  matrix  T  is  symmetric 
under  our  assumptions. 

For  small  deformations  X'.'e  have 

p  -  1  +  P   , 

where  P  is  small.   Then  up  to  terms  of  second  order 

pp "  5rf  I+P+P'"  ,    ip|  -t;.  1  +  [P]  ,    s^  -,-  [?  +  P'"] 

■f  '^(1  -  [P])(l  +  P  +  P'''')(2V^   +  2V^  ^  [P+P'""]  +  k^^I      (P+P""')) 

^1      ^1^1  ^2 

Z   2V^   +  2V^  (P+P'""-[P])  +  2V^  ^  [P+p"]  +  kM^    (P+P"")   , 
^1      ^1  ^1^1  ^2 

where  the  derivatives  of  V  =  V(s  ,Sp,s_)  are  taken  for  p  =  1  or 
s,  =  Sp  =  S-.  =  0.   If  we  require  that  the  stress  vanishes  for 
p  =  0,  we  have 

(3.8a)  V^  =  0   . 

^1 
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Put   for   s-    =  s^  =   S-,   =  0 

1  c  J) 

(3.8b)  i|V     g     ---  X     ,         i+Vg     =  tx     ; 

then  for  small  deformations 

(3.8c)  T  -i   a[P]  +  |j,(P  +  P'"')   . 

This  is  the  classical  linear  approximation  for  the  stresr.  strain 
relations  with  the  Lame  constants  \   and  [i.   The  constants  X   and 
[I   are  necessarily  positive  by  nature. 

For  later  purposes  it  vdll  be  important  to  have  the  ex- 
pansion of  U(p}  =  od+P)  in  the  neighborhood  of  ?  =  0  up  to 
terms  of  fourth  order  in  P. 

According  to  (3.7)  the  stress  matrix i^  is  expressible  in 
terms  of  the  symmetric  matrix 

Y  =  PP'"'   . 
I'Je  have 

(3.8d)  t     =    IyI"^^'^  2Z  2nV^  (y-1)""S' 

n=l    ""n 

where 

V  =  V(s^,S2,S3)   ,    s^  =  [(y-1)"]   . 

We  can  expand  x  into  a  Taylor  series  of  powers  of  y~1  ai^d  find 
from  (3.8a,b) 

Z    =    (1  -i[Y-l]  +...)(2V^      +   2V^    -    [y-1]    +  ij-V^    (y-1)    +    ...) 
'^  ^1  "1^1  ^2 


=  kr-1]  +  [lir-D  +  o((Y-i)^)     . 


2' 
This    agrees  vjith    (3. Be)    since 


Y-1   =   P  +  P'"'  +   Pp'""  =   P   +   P'""  +  O(P^)       . 


In  particular  for  y  =  1   and  arbitrary  dY 


'L  =0   ,     dr  =  ^[dY]  +  M-  dY 
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If  dr.  =  0  for  Y  =  1  and  some  dy  then 

0  =  |[dY]  +  i-i.  dY 

and  hence   also 

0  =    [0]    =    {l.\  +  ti)[dY]       . 

Since  X   and  p,  are  positive  this  would  imply  that  [dy]  =  0  and 

thus  also  dy  =  0. 

VJe  can  conclude  that  if  Y"1  is  sufficiently  small  then 

t  7^  0  for  Y~l  ^  0.      Otherwise  there  would  exist  a  sequence  of 

matrices  y  of  the  form  y  =  1+e  a  ,  with  scalars  e   tending  to 
'n  'n      n  n'  n       ° 

zero  and  matrices  CT  which  are  botmded  and  bounded  away  from 

n  "^ 

zero  such  that  the  corresponding  stress  matrices  "^  vanish.   Foi 

a  suitable  subsequence  the  o  vjould  conver.'te  tov.'ards  a  non- 

n 

vanishing  matrix  dy  for  which  dT  =  0,  in  contradiction  to  what 
has  been  proved. 


'The  matrix  y  =  pp "  has  the  same  invariants  as  the  matrix 
g  =  p"p.   Vanishing  of  y-1  means  that  p  is  orthogonal;  if  the 
strains  are  small  then  y-l  is  sr:all.   VJe  find  then  that  for 
sufficiently  small  strains,  i.e.  for  sufficiently  small  £-1,  the 
stress  matrix  X    can  only  vanish  for  orthogonal  p.   Since  by 
{2.1l\.)   T  and  q  vanish  simultaneously,  it  follov/s  then  also  that 
q  /^  0  for  sufficiently  small  strains  unless  p  is  orthogonal. 

lie  have  from  Taylor's  formula 

1     M  o^l   ^2  ^3 


U  =  V(s   Sp,s->)  =  2_^ I     -'.i'',  1  r  V.  .  .   s    Sp   s.- 

where 


V 


'S^=SyO 
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VJe  have   for  i   =  1,2,3 

s^  =    [(g-1)^]    -    [(P  +  P"   +  P"''>)^] 


il       r /T, .  T>"-\a, 


where 


<3-9)  V=^l^  [(P+P*)-(P*P)'']      . 


Then 


sjl   3^2   S33   =    (s^^+3^^)'l    (32,  +  s^^+s^2)    ^    (^3o^^21-*-^12-^^o3^    ^ 


^ik 


s.. 


1-    2-3-  ^ilki3aT^ 

where  the    siain  is   extended  over   all    sets   of  non-negative    integers 
^''lo'%l'''2o>''ll'%2'^3o'^21'^12'^o3^    ""^^^ 

(3.10)  a^^-^-a^^   =   i^,      a^^+a^^-Ha^^  =   i^,      ^3o''^21^^13'"%3   "  ^3* 
We   finally   arrive    then  at   the   foxTiiula 

a., 
s    ^^ 

(3.11)  u  =  rvi  i  ^      TT    ~^-"^i- 

^1^2^-3   i+k<3      ik* 

where  the  sum  is  extended  over  all  vectors  (a,   ..,,a„-)  whose 
nine  components  are  non-negative  integers  and  where  i^,ip,i-  are 
the  f'onctions  of  the  a,,,  defined  by  (3.10). 

By  (3.9)  the  expression  s  ^  is  a  form  in  the  elements  of 
the  matrix  P  of  degree  a+2p.   Hence  the  general  term  in  the  sura 

(3.11)  representing  U  is  a  form  of  degree 

(3.12)  N  =  JZ      (i+2k)a^^ 

i+k<3 

=  ^10  ■'2(a^i+a2^)  +3(a^^+a3^)  +^(^^2+0.^.^)  ^^a^^  +  Ga^^, 
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Each  set  of  non-negative  integers  satisfying  this  relation  fur- 
nishes a  homogeneous  term  of  degree  W  in  the  expansion  for  U. 
VJithout  attempting  to  determine  in  general  the  number  of  terms 
of  degree  N  it  is  easy  to  find  successively  all  solutions  a., 
of  (3.12)  for  i']  =  0,1,2,3,^'  ^^^   finds  in  this  way 

U  =  V     +  V-,   s,   +  V,   s  ,  +  i  Vo   s?   +  V  T^s^„ 
oco    loo  lo    loo  ol    2  2oo  lo    olo  2o 

1      1 

+  -r-  Vo    St    +  Vo    St   S  ,   +  V.,   St   S^    +  V  t   Stt   +  V  ^nS-j^ 

6   3oo  lo    2oo  lo  ol    llo  lo  2o    olo  11    ool  Jo 

+    TtT    V.  S7        +    -^    V-,         St"     S     T      +    TT    V^T      St      S^        +    Vtt^S-i^S,, 

2[).     i|.oo   lo        2     3oo   lo   ol        2     21o   lo   2o  llo   lo   11 

1  2  1  2 

+    Vt      tSt      S-,        +    Vtt      S     tS^        +    -^5    Vo        S     t     +    -:=-    V     o^S^ 

Icl   lo   3o  llo   ol   2o        2      2oo   ol        2     o2o   2o 

+  V   T    s    o  +  V      tSot    +    (terras    of  higher   order) 
olo   o2  ool   21         ^  ^ 

Without  restriction  of  generality  \-je   can  assume  that  V^„^  =  0. 

"        •-'  ooo 

Moreover  by  (3.8a,b) 

loo     S-,      '   ^  2oo      '   '  olo    ^ 
The  cubic  terras  contain  3  additional  constants 

6  3oo  *       llo  *       ool 
and  the  fourth  degree  terms,  [|.  further  constants 

^  =  -k   ^[.oo  '2=1  ^210  '   ^^^  =  V-^ol  '      ^^  =  \   Vo2o   • 
Vie   write 

(3.1)  u  =  u^(p)  +  u-'^(p)  +  iA(p)  +  ... 

where   U^'"(P)    is   a   form  of  degree   k.      Here 

U^  =  §   [P+P""]^  +  g   [(P+P")^] 

U-"^   =  ^    [P+P''''][P'*''P]    +  ^    [(P+P ■'■)?'■?]    +  A[P+P'"]^ 
+   B[P+P'"']  [(P+P'"')^]    +   G[(P+P"'')^] 
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U^  =  H  fP"P]^  +  77  [(P'^'P)^]  +  3A[P+P''']^[P""*P] 
o  ^ 

+  2B[P+p'''"]  [(?+?'"")?■'■■?]  +  3[P'"P]  [(P+P'")^]  +  D[P+P'"]^ 
+  E[P+P'''']2[(P+P")2j  +  P[.^+P'''][(P+P''')^]  +  G[(P+P"")^]^ 

+  3c[(p+p''')'^p'"p]  . 

For  any  matrix  ^  =  (4-i  )  ^®  ^^^   ^^®  symbol  D-  for  the 
differential  operator  defined  by 

(3.1^)  D^  =  ^,, 


Vie  have  for  independent  matrices  .^,yi  ,^,C,  .  .  .  and  for  p  =  1 

Plere  for  the  form  U^'CP)  of  degree  k  we  have  in  IT  (^  ,E,    ,  .  .  .  ,E,    ) 
^  polar  form  of  U^,  which  has  the  value  klU^(P)  for  E,^=£,^=.., 
=^  =P,  and  is  symraetric  In  t,    ,£,,..,  ,c,    .      It  is  simpler  to  write 
the  expressions  for  these  polar  forms,  v;hich  have  the  k-deriva- 
tives  of  U  i-jith  respect  to  all  the  p.,  as  coefficients,  than  to 
give  expressions  for  the  individual  coefficients.   We  find 
immediately  by  polarisation  from  the  expressions  derived  for 
the  U^(P)  that  for  p  =  1 

(3.15a)      D^D,U  =  X5^U,^)   =  T~  (r[^+€''""][^i  +  a''""]  +  ji[  (?+?"')  (v^+ h")  ]| 
(3.15b)      D   DD^U  =   U^(^,rj,^) 


J  Artro-tr  "■■!  r.,  "Vl    J.    Mil 


+    C[(^H-g'')(r,+v^")(^+^")]| 
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(3.l5c)   D^D  D^D^.U  =  u^(£;,r|,?,cr) 


T-rr 


;  \  ^  [E,'\][K''<y]   +  I  f^'^^'V] 


+  BU"''>i][(^+^'^'"")(^r  +  (T'"')] 

+  E  [C+?'''']  ['-,  +  <,■'•■]  [  (^+^'"")  (r  +(y""")  ] 

+  P[?+^''']  [(v,  +  '[■'■')  (?+^'"")  (0-  +  /)  ] 

+   3C[(^+4")(:+'/")^'V]l  . 

J 
Here  the  sums  on  the  right  are  to  be  extended  over  all  permuta- 
tions of  £,,f[    respectively  of  £,,r^,K   respectively  of  ^,n,^,cr. 
V,e   notice  in  particular  that 

^"^ikA,   .  ,   .f      ?;^U(p) 


^r.,u,'"^^'-  °  i^^iT^^,;.,'^'^''- "  '^"^" 


=  ^^ii'ikk  ■"   ^^^ik^iik  -"^iki^  • 

Having  the  expressions  for  the  derivatives  of  U  at  p  =  1  we  can 
easily  derive  expressions  for  the  derivatives  of  U(p)  at  any 
orthogonal  p  =  c.   VJe  have  indeed  by  {3 '2.) 

(3  17a)  r^'^ik^P^'l   £  v.   -  /"  ^^U(p)  "^   r   , 


^■p    1   ^ik'-rs   Ivp.,  ap  „/   ^ik 'rs 
•  p=c         "        -^p-c 

=  (D^D,^U(p))p^^ 
2 


=  6  ..  D  ...  c(pH   =  y^(c*?,o"'t() 
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2 

with  U  given  hj   (3.l5a).   Similarly  for  c  orthogonal 

(3.17b)       t^5D^D^U(p)]p^^  =  U^(c"''s,  c%,  c'-'^)   . 
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Chapter  IV 

The  Second  Variation  of  the  Strain  Energy.   Uniqueness 
at  Equilibrium  for  Prescribed  Boundary  Displacements. 


V/e  consider  an  equilibrium  position  x'  =  x   (x)  of  the 
material  maintained  by  certain  body  forces  P,  and  certain 
prescribed  tractions  or  displacements  on  the  boundary  S^  of  the 
region  R'  filled  by  the  material  in  the  strained  state.   Let 


ik'l) 


^l^.- 


;)x 


o' 


dx. 


'L  = 


?U(p) 


=>p 


ik 


/p=p° 


Then  by  f2.l5),  (2.16) 

.0 


(1^.2) 


3q 


ir 


^ 


|p°|Pi  =  0 


in  R 


^ir^r^S  =  X°^4°'  dS. 


for  X  on  S,  x 


o' 


on  S' 
o 


Let  x'  =  x' (x)  be  an  arbitrary  transformation  (not 
necessarily  corresponding  to  an  equilibrium  state).   We  compute 
the  difference  of  the  total  strain  energies  corresponding  to  the 
transformations  x'  and  x°  .   IVe  have,  expanding  U(p)  by  '-faylor's 
theorem. 


where 


'f 


U(p)dV  -  \\j 

■J  J  J 
R 


U(p°)dV  =  A  +  B 


B  =  ^ 


A  =  \\\(i^]  (p 


ik  -  Pik)  '^v 


■U(p) 


2/JJl5P.v^P 


R 


ik^  ^rs./  _^o 


(Pik-Pik^(P 


rs 


-p,3)dV 


p=p''+0(p-p  ) 


and  ©  =  0(x/  lies  between  0  and  1.  The  quantity  A  represents  the 
first  variation  of  the  total  strain  energy,  while  B  (at  least  for 
infinitesimal  x'-x   )  represents  the  second  variation. 


ko 


Using  the  equilibrium  conditions  {I4..2),    (1|..3)  we  have 


(i+.5) 


A  = 


Q,-i 


4i  -  4) 


ex. 


dV 


.0' 


JJjlp°|Fi(x.    -   x°")dV   + 


'?k«k<H-4')ds 


Assume  that  the  forces  acting  on  each  volume  element  in 
the  strained  state  are  prescribed,  i.e.  that 

P^  dV  =  P^lpldV 

is  known.   Let  moreover  on  a  portion  of  the  boundary  the 
stresses  acting  on  a  surface  element 


t.,  ^'dS.  =  q,,  £,  dS 


•ik^k 


ik'^k^ 


be  given,  whereas  on  the  remainder  of  the  boundary  the  displace- 
ments xt  are  Prescribed.   Vie  can  define  for  any  transformation 
x'  =  x' (x)  the  potential  energy  by 

(U  H'  '"/ 

p  =  !  !  u(p)dv  -  ijjp,|p|(x^ -x^)dv  -    j    qik^k^''i'''i^^^ 


1 1 


R 


where  p  is  the  Jacobian  matrix  determined  by  x'(x).   S,  is  that 
portion  of  S  for  which  the  stresses  are  prescribed  and  P  jp|dV 
and  q.,^  dS  have  the  prescribed  values.   If  on  the  remaining 
portion  of  S  the  displacem.ents  x'  and  x'   agree  we  have  for  the 
difference  in  the  potential  6ner3ies  corresponding  to  the 
transformations  x'  (x)  and  x'  (x) 


P-P-  = 


!  r 


R 


(U(p)-U(p"))dx  -  ill  Pjp|(x|-xt_°)dV  - 


R 


q.,  4,  (x!-x!  )dS 


=  B 


Since  B  is  quadratic  in  the  first  derivatives  of  x'-x'   we  have 
that  the  potential  energy  is  stationary  for  the  equilibrium 
state . 


See  Reissner  [51  for  a  different  formulation. 


kl 


The  second  variation  B  is  of  the  form 
(ij-.b)  B  =  ill  Q  dV 

with 


Jj 

R 


'^•^>    ^=H^iSfc)  <pi. -»?."■="» -p°-' 


rs  "^rs 


(i|..8)  p  =  p°  +  ©(p-p°) 


If  B  is  non-ne^^ative  for  all  displacements  x'-x°  of  a 
certain  type,  which  do  no  work,  then  the  total  strain  energy- 
will  be  a  minimum  in  the  equilibrium  state  for  all  those 
admitted  x' .   This  would  clearly  be  the  case  if  the  quadratic 
form  Q  in  the  Pj_jj-P^i,  is  positive  definite  or  semi-definite. 
Now  in  the  notation  (3.1ij.) 

ik.9)  Q  =  I  (D^Y"P"€=a=P-P° 

"p=p 

For  the   special  value  p  -  1  we   have   from   (3.l5a) 


=  AUJ^  +  ^  ((4  +  £.'')iC  +  ^)]  >  0 


where  B,   =  p-p°.   This  form  in  E,   is  positive  semi-definite.   It 
vanishes  whenever  ^  is  a  skew  symmetric  matrix.   There  is  no 
reason  to  assume  that  Q  will  still  be  semi-definite  when  p  is 
just  close  to  1  instead  of  being  equal  to  1,   In  spite  of  the 
expected  indefinite  character  of  Q,  there  are  cases  where  it  is 
possible  to  establish  that  at  least  B  is  non-negative.   This 
tjrpe  of  situation  has  been  studied  in  recent  years  by  Garding, 
Aronszajn  and  others  for  general  linear  elliptic  systems. 
Vie  observe  that  the  form 

(4.11)    xc^^c,,,,  +  t^(?ik«ik-^«iAi*  -  ^<5ik?ki-?ii?kk' 


kl 


The  second  variation  B  is  of  the  form 
(Ij-.b)  B  =  j  I  I  Q  dV 

with 


il^ 


(i|..8)  P  =  P°  +  ©(P-P°)   . 


o' 


If  B  is  non-negative  for  all  displacements  x'-x   of  a 
certain  type,  vjhich  do  no  work,  then  the  total  strain  energy 
will  be  a  minimum  in  the  equilibrium  state  for  all  those 
admitted  x'  .   'I'his  would  clearly  be  the  case  if  the  quadratic 
form  Q,  in  the  Pj^y'Piir  i^  positive  definite  or  semi-definite. 
Now  in  the  notation  (3.1i|-) 

(4.9)  «  =  I  <D.Y'P"€=n=P-P° 

p=p 

For   the    special  value   p  =   1   we   have   from    (3.l5a) 

(U.IO)       2Q  =    X[?]2   H-    y.[^^+^C]    =   ^^ii^i^ir    +   l^(^ik^ik-*"^ki^ki^ 

where  4  =  p-p°.   This  form  in  E,   is  positive  semi-definite.   It 
vanishes  whenever  ^  is  a  skew  syminetric  matrix.   There  is  no 
reason  to  assume  that  Q,  will  still  be  semi-definite  when  p  is 
just  close  to  1  instead  of  being  equal  to  1.   In  spite  of  the 
expected  indefinite  character  of  Q,  there  are  cases  where  it  is 
possible  to  establish  that  at  least  B  is  non-negative.   This 
type  of  situation  has  been  studied  in  recent  years  by  Garding, 
Aronszajn  and  others  for  general  linear  elliptic  systems. 
We  observe  that  the  form 

(It. 11)  xc,,?,^,^  +  ^(?i^5ij,^-q^?„i)  -  ^(?ik«kl-5ll?kk) 

=  (XH.^x)(C,,)2.^?,^q^ 


is  positive  definite  in  Ey    since  X  and  ^   are  positive  constants. 
The  same  holds  then  for  all  quadratic  forms  with  neighboring 
coefficients.  Assuming  U  =  U(p)  to  have  continuous  second  de- 
rivatives with  respect  to  the  p.,  there  exists  then  an  e  only 
depending  on  the  choice  of  strain  energy  function  U  such  that 

=  ^D,D^U(p)  -  ^(^k^ki-^li^kk)  >  ° 

f  7^  0,  (p-"-l)(p-l)  <  €  . 

Hence 

ik.iz)  0  >  ^  ((Pi^-P?^)(P^l-PM)-'Pll-P?l'(Pkk-PkkV 

=  I  ?%  (<pik-p?k'<4-£')-(Pkk-pSk)(^;-4'>) 

for  p  7^  p°  and 

(l|..13)    [(p-^-l)(p-l)]  <  e  ,  [(p°''''-l)(p°-l)]  <  e  . 
It  follows  that  under  the  assumptions  (lt..l3) 

E  ^  t-;^  ('pik-p^kX^^i-^'S' )  -  <Pkk-Pkk"='i-4' > ) ^1  ^^ 

unless 

Pik  =  Pik   ^"  ^  • 

t      o ' 

If  the  displacements  x  and  x   are  equal  on  the  boundary 

and  if  p  and  p  lie  in  the  e -neighborhood  of  the  identity  de- 
scribed by  (ii..l3),  we  find  that  B  >  0  unless  p=p°  in  R.   The 
identity  of  p  and  p°  implies  that  the  xl-x?  are  constant  and 
hence  vanish,  since  they  vanish  on  the  boundary  S.   This  yields 
the  theorem: 

In  the  absence  of  external  volume  forces  any  change  from 
the  equilibrium  state  that  leaves  the  boundary  points  fixed 
leads  to  an  increase  in  total  strain  energy,  as  long  as  only 
states  in  an  e -neighborhood  of  the  identity  are  considered. 
Here  the  e -neighborhood  of  a  transformation  with  Jacobian  matrix 
p°  is  defined  by 

{(?%"■")  (P-P°)]  <  e. 


>w 


•  i-     y}r'n 


J  <:  Si-. 


<  V  ; ;  i  i  "J. 


e.-'j-   ■. 


f-ofp    <*•    ""t  f     !*'" p 


k3 

Since  then  for  given  boundary  displacements  any  state  dif- 
fering from  a  given  equllibri'jm  state  (in  an  e-neighborhood  of 
the  identity)  corresponds  to  a  larger  total  strain  energy,  we 
arrive  at  a  uniqueness  theorem: 

For  given  displacements  at  the  boundary  there  can  exist  at 
most  one  equilibrium  state  in  the  e-neighborhood  of  the  identi- 
ty. Here  e  depends  exclusively  on  the  choice  of  strain  energy 
function. 

It  would  be  desirable  to  havie  a  uniqueness  theorem  apply- 
ing to  transformations  x'  not  restricted  to  the  e-neighborhood 
of  the  identity.   It  is  clear  that  the  validity  of  such  a  theo- 
rem would  depend  on  the  behavior  of  U(p)  for  matrices  p  ranoved 
from  the  identity. 


Chapter  V 
Two  Dimensional  Bending  of  a  Thin  Plate 

In  this  chapter  we  only  consider  plane  deformations  of 
the  form 

\  ^  m  X  )  X-.  —  X-i  \  X"!  J  Xp  )      f     Xp  —  Xp  V  X-j  Xp  /  f  X^  —  X^  . 

Then  the  Jacobian  matrix  p  is  of  the  form 
(5.2) 


) 

/Pll 

Pl2 

p  = 

\:" 

P22 

0    i 
1 

0 

1/ 

It  follows  for  the  matrix  y  =  PP  that 

^13  "^  "^23  "^  "^31  "  *^32  ^  ^'  ^33  ^  '^' 

Since  by  (3.8)  the  matrix  ,"'  is  a  function  of  y  we  have  then 
also  .  ^  ■, 

^13  "  ^23  "  "31  "  "32  =  °'  -33  "  2|prVg^(s^, 82,33). 

Thus  for  q  =  |p  I  "^^  (p '"■)"■'■ 


(5.3) 


"^13  "  "^23  "  ^31  "  "^32  =  °'  ^33  ""  ^^s^^^l'^2'^3^ ' 


i  X 


For  a  matrix  p  of  the  special  form  (5.2)  the  s.  are  determined 
completely  by  Pi  t  »P-i2»P2i»P22*   ^^^^e  those  quantities  are  known 
we  can  find  all  of  the  'K,  and  q. ,  explicitly. 

Using  (5.3)  and  the  fact  that  all  Pii^-j^^k  ^^®  functions 
of  x,,X2  only,  the  equilibrium  equations  (2.l8)  in  the  absence 
of  volume  forces  reduce  to  the  two  conditions 

i^.k)  ~^-  +  —^  =0   for  i  =  1,2 

-■^1    "^2 
In  addition  we  have  the  compatibility  conditions 

(5,S)  — — ^  =0  for  i  =  1,2. 

The  matrices  p  and  q  are  connected  by  the  relations 

IK    r-Pik 
In  forming  the  q.,  with  i,k  =  1,2  the  matrix  p  can  already  be 
assumed  to  be  of  the  special  form  (5.2), 

We  consider  now  the  cas3  of  a  medium  which  in  the  unstrained 
state  fills  the  region  with  rectangular  cross  section 

(5.7)  0  <  x^  <  a  ,  -h  <  X2  <  +h  . 

J, 

The  boundary  requirements  on  the  "horizontal"  sides  x^  ~   ~^ 
shall  be  the  absence  of  tractions,  or  by  (2,21) 

(5.8)  q^2  ~  ^22  ~  ^   ^°''^'  ^2  ~  "^  ' 

Boundary  conditions  on  the  vertical  sides  x,  =  0,  a  would  have 
to  be  imposed  as  well  but  will  not  be  specified  at  present. 

One  trivial  family  of  solutions  is  obtained  by  taking  for 
X,  ,X2  linear  functions  of  x.,,X2.   The  matrices  p  and  q  will 
then  have  constant  elements,  so  that  (5»U)»(5.5)  are  satisfied 
automatically.   The  constant  matrix  p  has  only  to  be  chosen  in 
accordance  with  the  boundary  condition  (5»8)  that 

q^2(p)  =  q22^P^  ~  ^' 

¥e   are  interested  in  other  types  of  solutions  correspond- 
ing to  a  buckled  state  of  the  medium  for  small  thickness  h. 


v.;.f*- 


. '.: 


>iL 


U.^'      ' 


..      SSo  £ 


We  assume  that  there  exists  a  solution  of  (5«l<-)  »(5«5) » (5«6) , 
(5»8)  which  depends  analytically  on  Xp  and  h.   The  matrix  p  can 
then  be  expanded  into  a  power  series 

(5.9)     p  =  2_  ■xhrv''^   x^hP 

where  ,  ^ 

ap    ap,   ,    f  ^''"'^P(Xi,X3.h) 
p  ^  =  p  ^'(x,  )  =  I  -— ^ 

Because  p  is  of  the  form  (5.2) 

(5.10a)   p°|  =  p^j^  =  0  for  a  +  p  >  0;  i,k  =  1,2,3. 

Similarly 

(5.11)     q  =  q(x^,X2,h)  =  ^     _^^  q^P  x^  h^ 

where  the  matrices  q  ^  depend  on  x,  alone.  Here 

The  differential  equations  (5.i4-) ,  (5.5)  become 

,^  ^>    '^  ^11  _^   a+1  p   ^  ^  Pi2    a+1  p   ^   „   .    ,  ^ 
^^•^2)    -diT"  ^  ^12    -  0,  -^3^  -  p^^  P  =  0  for  1  =  1,2. 

The  boundary  conditions  (5.6)  yield  the  additional  equa- 
tions ap  ap 
(5.12a)    ^  ^=  Z:  (-1)''^=0 

for  i  =  1,2,  and  y  =  0,1,2,... 
Adding  and  subtracting  these  two  relations  and  expressing  q.p 
in  terms  of  q?r  by  (5.12)  we  arrive  at  the  two  equivalent  con- 
ditions 
(5.13a)    (J)q°^Y  ,  (pq22Y-2  ^  (Y)q^2Y-U  +...  =  q 


(5.13b)    (^;^)q?i^  ■»-  (^3^)qii^"^  +('^;^)q^lY-^  +...  =  const.=  C 
for  i  =  1,2   and  y  =  0,1,2,... 


n  '■'■ 


1^6 

Equations  ( 5.12 ),( 5.13a, b )  are  supplemented  by  the  re- 
cursion formulae  obtained  from  (5.6)  that  connect  the  q 
with  the  p  '^^.  We  have  first  of  all 


(5.lUa) 

(5.1i|b) 

Xp=  h  =   0 

_,'  '1  q       -^rs  "^ 

l^Prs  ■'"'Z    '      ..  i^Prs/   00  ^^ 

X2=h=o  *■     P^^P 

=  '°  10  ^'  .  oo   • 

p    p=p 

Similarly 

(5.11+c)    q°l  =  (D  ^^  q)    ^^ 

P      P=P 

(5.1l|d)    q^°  =  (D  i^D  3^^  q+  D  20  q)    oo 

P    P        P      P=P 

(5.1i^e)    qll  =  (D  ^^  D  ^^  q  f  D  3_^  q)     oo 

P    P        P 

(5.1i;f )    q°2  =  (D  ^1  D  ^^  q  +  D  ^2  ^) 

P    P        P      P=P 
etc. 

For  the  quantities  of  lowest  order  we  have  from  (5.13,a,b) 
for  Y  =  0 

(5.15)  qll   ^  ^'  *^il  "  const.  =  C°  for  i  =  1,2 

The  p*^*^  would  then  have  to  be  determined  from  relations  (5«ll4-a) 
For  kno^^)n  p?°  q.?  the  values  of  the  q.p  and  p.?  follow  from 

(5«12)  for  a=p=0.   The  remaining  components  p.p  would  follow 

from  the  relations  (S.ll+b)  from  the  known  q^p*   Prom  (5.13a,b) 
for  Y  =  1  we  have 

(5.16)  q?p  =  0,   q.,  =  const.  =  p  ^i  ^^^  ^   ~   •'•»^* 


■.,v  -      p,  ■' 


For  known  p°°,q°  we  have  In  (5.lUc)  a  system  of  linear  equa- 
tions for  the  Pj^i..   If  this  system  has  a  unique  solution  we  can 
proceed  in  the  same  manner  and  get  successively  all  p   ^ ,    q   ^ 
If  p  °  has  been  determined  and  the  constants  cT  are  known. 

The  values  of  the  constants  cT  are  determined  if  suitable 
boundary  conditions  are  prescribed  on  the  vertical  sides  x,=0,a. 
These  constants  are  indeed  determined  uniquely  by  the  resultant 
stress  force  applied  to  the  vertical  faces.  We  have  from  (5.13b) 

(5.17)   ■:  qii(x^,X2)dx2  =  ;  JZ     ahr  ^ti^^i^^P'^  <^^z 
-h  -h  ^'P=° 


l^_^  rmiW  <^x,)(i.(-i)")h-P^i 


For  x-,=  a  we  have  by  (2.16)  for  the  components  of  the  result- 
ant of  the  stresses  acting  on  the  face  x,=  a  (per  unit  of  length 

of  X-.) 

-^      "  h 


(5.18)      Sr^i^'^l   ds'  =  ,^  q.-^(a,X2)dX2 

-h 
If  this  resultant  stress  is  given  for  each  h  the  constants  C. 
are  determined. 

The  procedure  outlined  determines  for  known  p°°  and  cT 

uniquely  all  p   ,  q   and  hence  also  p  and  q,  provided  equation 
(5.11+c)  can  be  solved  for  p°  ,  equation  (S.li+f)  for  p°  ,  etc. 
This  is  the  case  as  long  as 

^^F.    ^)    oo  =  (  ^]  ^rs 

p=p     V  ^rs  ^p=pOo 

cannot  vanish  for  f  7^  0.   The  interesting  cases  corresponding 
to  buckling  occur,  when  this  condition  is  violated.   These  are 
the  cases  where  the  determinant  of  the 

?^ik  _  ?/v 


'^Pik   ^Pik^Prs 
vanishes  for  p  =  p°°  , 


^Vf    ^-t 


■^-1 «' 


C.-''V 


".9'-£.1 


._,  V    .''>  ';!• 


t:  :  ■•      •!.■ 
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This  exceptional  situation  always  arises,  when 

(5.19)     q°°  =  0, 

i.e.  when  the  stresses  vanish  to  lowest  order  in  h.   Condition 
(5«19)  has  been  shown  in  Chapter  3  to  be  equivalent  to  ortho- 
gonality of  the  matrix  p°*^. 
Let  now 

(5.20)     p°°  =  c     where  cc'"'  =  1. 
By  (3. 17a), (3. 15a) 

p=c  ^ 

p=c 

=  \J^{o"\,c''rj)   =  \[c''e][c"V]+ti[c'"^c'''/r +?'^;]. 
The  relation 

(5.22)  (D^q)p^^  =  0 

for  a  certain  f  is  equivalent  with 

U^(c'"'^,c'"'ir})  =  0  for  all"";. 
For  any  matrix  a  the  expression 

^^  -^  =  ^kl^lk 
vanishes  identically  in  ^j   if  and  only  if  a  =  0.   Hence  (5.22) 
is  equivalent  to 

\[c''V]c'^+  |a(c'"?c'V"')  =  0 
or  to 

(5.23)  Mc\]    +  ^{c'\-^E,%)   =  0. 

Taking  the   trace   of  relation  (5.23)   we   find 

(3A  +   2^)ic'\]=   0 
and  hence    (for  positive   \,\i} 

lc'^^]=   0 
and   then 

{c'\)   +   (c'-V)"  =   0. 


(  V 


(  ■..-■; 
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Thus  (5.23)  la  equivalent  to  skew  symmetry  of  c''^,  i.e.  to  the 
existence  of  a  matrix  f\   such  that 

(5.2U)   ?  =  cA.  ;  rt-=  -a 

Whenever  (5. Si;)  holds  then  for  any 

(5.2L|.a)    [(D^q)p^^T|'^"]  =  0. 
Assume  then  that 
q°°=  0 
which  is  consistent  with  (,^.15).   The  matrix  p°°  is  then  orth- 
ogonal.  Since  p    is  of  the  form  (5»2)  also  we  must  have  an 

angle  6  =  6(x, )  such  that 

;cos  6     -sin  6     0  \ 

(5.26)  p°°  =  c  =     !  sin  e      cos  6     0  i 

\    0  0     1/ 

Ve  make  the  further  important  assumption  that  6(x, )  is  not  con- 
stant, thus  excluding  trivial  solutions.  ,  i 

Along  the  line  Xp  =  const.  =  0  we  have  tan  6  =  . — -  so  that 

'^  dx| 

e  gives  the  inclination  of  the  "middle  surface"  in  the  strained 

state  in  the  limit  for  h  =  0.   It  is  our  object  to  derive  a 

differential  equation  for  6. 

The  q?.  for  i,k=l,2  have  constant  values  given  by  (5.16). 

By  (5.1l;c),(5.21) 

(5.27)  q^j^  E^^   -  ,  -^—   r^^p^^  -  U  (c  p   ,c  ?)  -  0 

when  r    is  of  the  form  (5.2L|.).   Take  for  CX   the  matrix 
/  0     -1     o\ 

(5.27a)  I .:  =  I  1      0     C 

0      0     0/ 


\ 


so  that        \ 

/'-sin  e  -cos  e  o\ 

i  \         H. 

(5.28)  r   =  c  n «j   cos  6  -sin  6   0  i  =  f| 

\      0        0   0/ 


i>o?-': 


h'.'  \.'0 


■<  ♦ 


i     V    .        I 
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Then  by  (5.27) ,(5.16) 

for  all  X,.   Since  9  is  not  constant  it  follows  that 

ci  =  ci  =  0 

and  hence  (see  (5.11a)) 

(5.30)    q°j^  =  0  except  possibly  for  i=k=3. 

There  is  no  need  to  determine  p°  completely.   ^«'e  have  from 
(5.30),(5.iOa),(5.il;c),(5.2i) 

0  =  <l   P?k  =  (D  ^,D  ,itJ(p))    =  u2(cV\c%^l). 

P    P 
By  (I+.IO)  the  form  in  ? 

U^(^,^)  =  A[^]2  +  |[(H  +  5'')(£*+  r)] 

is  semi-definite  and  vanishes  only  for  skew  symmetric  E*      Hence 

(5.31)     (cV^)  +  ic%'''^r  =  0. 

For  q°  we  have  from  (5.1^f)  and  any  £ 

P   P      P    P=c 
For  £  given  by  (5.28)  it  follows  from  (5.2i4.a)  that 

[(D  pq)   C''i  =  KD  q)   p^^'^]  =  0. 

p    p=c         ^  p=c 

Consequently,  using  (3.17b), 

^^•32)     q?k^ik=  tq°V]  =  [(D^;,D  ^,q)    ^1 

P   P    P=c 

=  ^°  ol°  oA^(P))   =  u3(cV^.cV^,c*4) 

p   p   '^     p=c 

_  tt3/  '^  ol   -"■  Ol  .-,  V 
=  U-^(c  p   ,c  p   ,il). 

The  general  expression  for  IJ-^  has  been  computed  in  (3.15b). 
Since  here  all  three  arguments  c  p   ,  c  p  ,  sl    are  skew  sym- 
metric  and  each  term  in  V^   contains  at  least  one  factor  depend- 
ing on  the  symmetric  part  of  one  of  the  arguments,  we  cleaiiyhave 


o\ 


f.-v.<i 


•  L-^-   ■' 
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'J     Cri.^ 
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(5.33)         q?^  ^1,  =  u3(cV\cV\n)  =  0. 

By    (5.13a,b)    with  y=2 

(533a)  q°|  =    -q2°    ,    q°^   =    -  i(  qf °    -  C^)      for   1=1,2 

Hence  (5.33)  can  be  given  the  form 

=  ^(C^  £^^  +  C2  ?2l)- 
'"'e  can  write  this  relation  in  the  form 

Here  by  (5.12)  for  a=l,p=0 

20  .   _   '^  <^11  , 
^12  ^12      dx^  ^12 

whereas  by  (5.1tj.d)  (in  analogy  to  (5.32)) 

4i  ^k  =  (D 10°  loV^p))  , 

p   p   ^     p=c 


=  U-'  (  c  p   ,  c  p   ,-■■). 


Hence  our  relation  {S»3h)   goes  over  into  a  relation  only  involv- 
ing the  matrices  p  °,q  °: 

dq^° 
(5.35)    C^  q^  =  -2  -^  ^2  ^  U^(cV°,  cV°,i"). 

■5 
In  the  expression  (3.l5b)  for  U-^  the  terms  containing  the 

higher  elastic  constants  A,B,C  make  no  contribution  here,  since 

P.  is  skew  symmetric.   Hence 

u3(cV°,cV°.:'  )   =  2X[cV°][^''cV°] 

+  2ti[cV°"  c)/VcV°] 
=   2X[c^^°][?Vl    ^2ix   [c^^p^%^°"c)5V°] 
=   2u2(cV°,?V°) 

=   2(D   i^D  ,^   U(p))  =   2[qlV°'"'^c-=^-] 

P         cC  p  P=c 


•      OM* 


S2. 
For  the  matrices  ^c, a"!  given  by  (5.28) ,{ 5.26) ,  (5.27a)  we  have 

^c""  =  O 

Hence      o  -  t  -^  ..  i  t   i  ''^ 

U3(c'-pl°,c"pl°,;))  =  2[q^°pl°  r-] 

-  2(^lk  P2k  -^2k  Plk^ 
The  differential  equations  (5.12)  yield  for  a=p=0 

nlo  -   -^1  -  o    lo  -  ^12  _  lfi2 
%2  ~  "  dx^   ^'   Pll    dxj_  "  dx^ 

Hence 

U  (c  p   ,c  p   ,.(;  )  -  2(q^^  -^~  -  ^^^     -^) 

_  p,  lo  ^^11  ,   ol  ^^21  ,  _  p  ,lo  '^'^il 

On  the  other  hand 

P  ^^-11  ,        _  p  <^^il  , 
"  ^  dx~  ^12  ~  "^   ~dy~  °il  * 

Consequently  relation  (5*35)  takes  the  form 

<5-36)    C?  ?ll  =  2  3%  <<JU  =ll' 

It  remains  to  evaluate  the  term  on  the  right  hand  side  of 
(5.36)  in  terms  of  expressions  involving  6.   i''e  have  for  any 
matrix  '') 

[q   ;  ]  =  u^(c  p  ,c  ;) 

=  X[c  p-'^lLc  'J  +ti  [(c  p-^"+  p-""  c)c  -'/]. 
Since  T(  is  arbltra]:^  this  implies  that 

ql°  =  X[cV°]  c+tx  (p^°+  cp^°"c). 
Multiplying  this  relation  from  the  left  with  c  we  obtain 

*  lO    -V  r   """  lOi   .    /   '"*  lo.    lo    \ 

c  q   =  A[c  p   ]  +  M-(c  p  +  p   c) 
or  component  wise 
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53 
For  i=k=l   and   i=k=2  we   obtain   respectively 

°rl  ^r?  =  McV°]   ^2^  0^^  pJ;J 
0  =   \[cV°]    +2ixc^2  Pr2 

=  X[cV°l   +2,0^3   pl°   -2^0^^  j>ll   -   2^c^3  pj;3° 

=  (X+2tx)[cV°]   -2m.c^^  pJ;J    . 

Eliminating   [c"p   °]   bet^iieen  the   two   equations   we   find 

^rl  ^rl  \+2ti        °rl  ^rl 

^  Uu(X4tjt)  ^  Pr2 

X+2|x  rl     dx. 


\+2'iJi        '^rl  ^r2  dx-^  X+^M-         dXj_    ' 


Substituting  this   expression  into   (5.28)   we   obtain  for 
6  =   6(x,  )    the   differential   equation 

(5.37)   ^^;M  dfe  ,^2  3.^  e  +  cf  cos  ©  =  o. 

The  coefficient 

(5.36)       D=|^^^h3 

is  the  "flexural  rigidity"  of  the  plate.   (See  Love  [  1]  ,p.[|.6U) . 
The  constants  C.  xrjere  defined  by  (5 •33a) 


'1 
'i  "  ^11  •  ^^11 


-.2    2o  .  -,  o2    „    .  ^  o 
C.  =  q. T  +  3q.T    for  i=x,2, 


The  terminal  force  (per  unit  of  length  of  x^)  acting  on  the 
face  x^  =  a  is  by  (5.lS) , ( 5.17) 


\   'ik  -^k  '^^      tr=n   ( Y+1 )  l   i 
x,=a 


This  is  also  the  total  force  acting  on  any  cross  section 

x,=  const,  in  the  strained  state.   Since  in  the  case  considered 


»■:• 


J  J 


•  ■■i,-  I 


;.  } 


>Vi"    -j-xj    -^l 


Here  the  C.  and  C.  vanish,  we  have  for  the  components  of  the 
terminal  force 

^  cj  h^  +  O(h^),   (i  =  1,2). 

The  component  of  this  force  vertical  to  the  middle  surface  in 
the  strained  state  (the  Shear  force  in  a  cross  section  x,  = 
const)  is  then  given,  up  to  terms  of  higher  order  In  h,  by 

T  =  l-(-C^  sin  e  +  C2  cos  6)  h^  +  0(h>). 

Hence  (5«37)  yields  up  to  terms  of  higher  order. 

d^e  _   T 
ax. 

Since  p°°  Is  orthogonal,  we  have  also  up  to  terms  of  higher 
order  In  h  for  Xp  =  0 

ds'^  =  dx^^  +  dx2^  =  (p^^^  +  V2-^)   ^^i^ 

=  dx^^d  +  0(h)) 

Hence  up  to  terms  of  higher  order 

D  S_^p  +  T  =  0. 
ds  '^ 

By  a  suitable  rigid  motion  we  cm  always  bring  about  that 
the  terminal  force  acting  on  the  end  x,  =  a  is  parallel  to  the 
x,-axis  and  has  components  (-R,0),  where  R  is  positive  in  case 
of  compression.   Then 

J  C^  h^  =  -R,  c|  =  0 

so  that  the  differential  equation  for  6  takes  the  familiar  form 

D  ™-  +  R  sin  e  =  0. 
dx"^ 


(vr,  O  &■: 
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Chapter  VI 
Buckling  of  thin  plates. 

We  consider  a  solid,  which  in  the  unstrained  state  is 
restricted  to  a  portion  of  the  parallel  slab 

-  h  <  x^  <  +  h. 

The  faces  in  the  strained  state  corresponding  to  x^  =  1  h  shall 
be  free  of  tractions,  i.e.  we  assume  the  boundary  conditions 

(6.1)  q^^  =  0   for  x^  =  t  h. 

Let  the  solid  be  in  equilibrium  under  the  influence  of  certain 
lateral  forces  with  external  volume  forces  absent.   To  derive 
differential  equation  for  the  "middle  surface"  for  h  — ^  0,  we 
proceed  exactly  as  in  the  preceding  chapter.  We  assume  that 
there  is  a  family  of  solutions  depending  analytically  on  x, 
and  h,  for  which  p  and  q  are  given  by  power  series  expansions 

?CD     X,   h"^   „n. 

(6.2)  P  =  ZZ   -|t0T-  P°-^i^,,^2^ 

a  ,  3 
OP   X,  h^   _R 

(6.3)  q  =  EI  -traT-     Q^^Cx^.Xg) 

The  differential  equations  of  equilibrium  (2.18)  become 

(6.1+)     qfl^^   +  qjP^2  "^  ^i3^  ^  =   ^        ^°^   i=l,2,3 

//    rr\  a3  a6  a3  a+1   3  aS  a+1   3 

(6.5)  p^g^i  =  p^l^2    '   Pi3,l  =  Pll  '   Pi3,2     =  Pi2 

for  i=l,2,3.      In  addition  we  have   from  the  boundary   conditions 
(6.1)    for  Y  =   0,1,2, .. . 


a3  ap 

,       ^%1  .0  5—     (-1)^     -^  =  0 

a+3=Y    ^^*  a+3=Y                   ^•'^• 

in  analogy  to  (5.12a).  The  result  of  adding  these  two  relations 
can  be  written 


.*-' 


ranrf^ 


V.'M' 


r.^; 
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(6.6)    (V,,o^Y  ^  (r„2^r-2  ,  (r„U^r-4  ,  ...o. 

If  we  subtract  the  two  relations  and  use  (6.i;)  we  obtain 

To  lowest  order  we  have  from  (6. 6),  (6. 5)  the  nine  equatloia 
(6.8)    ,-  =  0,  ,-_^  .  ,00^^  ,  ,^   poo^^  .  pOo_^  ,  „ 

for  i=l,2,3 

where  the  matrices  q°°  and  p°°  are  connected  by 

(6.9)     q°°  =  q(p°°). 

Ve  choose  again  among  the  solutions  of  the  system  (6.8), 

(6.9)  the  ones  corresponding  to  vanishing  stresses 

(6.10)  q°°  =0. 

Then,  as  we  know,  p   must  be  an  orthogonal  matrix.  We  can 
easily  find  the  most  general  p°°  =  p°°(x,,Xp)  which  is  orth- 
ogonal and  also  satisfies  the  last  group  of  relations  in  (6,8): 
The  equations 

Pl2,l  -  Pil,2  =  °'  i=l»2,3 
are  equivalent  to  the  existence  of  three  functions 

w,(x,  ,X5)  =  x.'°°  such  that 


'l\^l*^2- 


^Ik  ~  "l  k  ^°^  i=l,2,3  and  k=l,2. 


Orthogonality  of  p°°  Implies  that 

Pji  Pjk  =  ^ik       ^°^  i,k=l,2,3. 
For  i,k  =  1,2   this  yields 

Wj  j^  p°°  =   0       for  k=l,2 

This  means  that  the  vectors  (w,,W2>w-i)  =  (x|°°,X2°*^,xl°°) 

for  varying  x^,X2  describe  a  developable  surface  referred  to 
the  isometric  parameters  x, jXp.   Since,  moreover 


■°°  =  0   for  k=l,2 


we  see  that  the  quantities  p?.^  are  the  components  of  the  unit 
vector  normal  to  this  developable  surface. 

We  shall  only  consider  here  the  case  where  this  limiting 
middle  surface  is  the  plane  xi  =  0,  and  where  in  fact  the  whole 
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transformation  is  the  identity  to  lowest  order,  i.e.  we  take 

(6.11)  P~=6,^ 

This  is  consistent  with  all  relations  (6.8), 

ve  proceed  to  the  determination  of  p°  and  q 
We  have  from  (6. 5), (6. 6), (6. 7) 

(6.12)  ^ol  _  n  „ol   .  ^ol   _  ^   ol     „ol   _  f^ 

for  1=1,2,3. 
Moreover,   by   (5.1[|.c) 

q°^  =  (D  „,  q) 

P  P  =  l 

and  hence  for  an  arbitrary  matrix  E, 

^Ik  ^ik  =  f^°^  ^"'^  =  ^^  ol  ^^  ^^P^)  _  =  U^(P°^^) 

Since  this  expression  does  not  change  if  E    is  replaced  by  E 
the  matrix  q°  must  be  symmetric: 
(6.11;)     q°J  =  q°^  for  i,k  =  1,2,3 

We  then  have  from  (6.12) 

(6.15)  ql]   =  q°^  =  0   for  1=1,2,3 

/A  ■n'\  ^ol         .    ^ol   _  ^   „ol    .^ol   _  ^   ^ol  _  „ol 

(6.16)  q^^^^^  +  q3^2,2  "  °'  ^21,1  +^22,2  "  °'  ^12  "  ^21  " 

The  general  solution  of  equations  (6.16)  is  of  the  form 

(6.17)  qll   =  4^22  '  ^12  =  ^21  =  -  ^,12  '  ^22  =  ^11 
with  a  suitable  function  ^(x^jXp).   (Airy  function). 

We  can  show  that  4  has  to  satisfy  the  bi -harmonic  equation. 
For  this  purpose  we  write  Identity  (6.13)  component  wise 

(6.18)  q°J  =  X[p°M  6^j^  +  p^ivll   +  p°J). 
Ve  observe  that  by  (6.17) 

//  ,o\    .  2 1    „ol      „ol      „ol     .  _ol 

(6.19)  /^.  4  =  qii^22  -  ^12,12  "  ^21,21  ^   <^22,11 
where  -.^   denotes  the  Laplace  operator 

.  2     . 2 

/x.  =  ^c_^  +  _^ 

:5X,       -  X2 
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ve  have,  quite  generally,  if  f  and  2^  are  matrices  connect- 
ed by 

(6.20)  ^  =  (D^  q(p)) 

^      p=l 
or 

^ik  =  ^^rr  ^k  '   ^^(^ik  ^  ^ki^ 
that 

(6.21)  q^^  =  K^^     for  all  l,k  =  1,2,3 
and  that 

(6.22)  ^11,22"  ^12,12"  ^21,a  "^^22,11"  31+2]!^  ^  ^ll'^^22"^^33^ 

=  2ti(f^^^2-  ^12,1^2  "  ^^^^^21, 2  '  ^22,lh 

Applying  this  identity  for  p   =  p°^ ,    ^  =  q^-'",  we  find  from  (6.17), 
(6. 15), (6. 12)  that 

^1  -3x^1:)^^^-'^  =  ' 

and  hence,  since  X  and  \i   are  positive, 

(6.23)  l)?^  =   0. 

In  addition,  we  have  from  (6.15) , ( 6.l8) ,( 6.12)  that 

r^ol   -  ^ol     ^ol  .  ^ol  _  p,  ^ol  .   ol  _  « 
^31,2  -  P32,l  '  Pl3  ^  P3I  ~  °'  P23  ""  P32  -  °- 

Hence  there  exists  a  function  v  =  v(x,,X2)  such  that 
(6.2[^)    p^^  -  p^3  -  v^^,  P32  -  -  P23  -  v^2 

The  scalar  functions  4  and  v  essentially  describe  the  matrices 
p*^  ,q°  .   ^"'e  have  already  the  bi -harmonic  equation  for  (|>.   A 
differential  equation  for  v  follows  from  consideration  of  the 
terms  of  next  order. 

^'^e  have  from  (5.11^b) 

q^°=  (D     q(p)) 

P  P-1 

or 

(6.25a)   q^°  =  X  p^°  d.^   +  ti(pjg  +  p^°)   for  i,k  =  1,2,3. 

Moreover,  from  (6.U),(6.5)  for  a=p=0,  using  q°°  =  0,  p°°  =  1: 

(6.25b)    q\°   =  vll   =  Pi2  =  0  ^°^  ^  =  1,2,3. 

It  follows  immediately  from  (6.25a,b)  that  all  elements 


D/  ;  i 


.^,/ 


*  ^'■■ 


•>  n 


59 

It  follows  immediately  from  (6.25a,b)  that  all  elements 

lo    lo    .  , 
Pik  »  ^ik  va^is^ 

(6.26)  p^°  =  q^°  =  0. 
Then  by  (5.lUd) 

P     P=l. 
Since  apaln  from  (S.l;) ,  ( 6.5)  for  a=l,  (3=0  using  (6.26) 

qi3  =  PiJ  =  Pi2  =  0  foJ^  i=l,2,3 
it  follows  in  the  same  way  that 

(6.27)  q2°  =  p2°  =  0. 

Generally  the  assumptions  q°°=  0,  p°°=  1  imply  that 

(6.27a)    q^°  =  p^^^  =  o  for  a=l,2,3... 

Equations    (6. 6), (6. 7)    for  y=2  then  yield 

(6.28)  q°^  =   0,    q°^^^   +   q°|^2  =   °      ^°^  1=1,2,3. 
Moreover,    by   (5.1i|f) 

(6.29)  q°2  =  (D  ^-^  D  ^;L  ^  +  D  ^2  q)     . 

P    P        p     p=l 

Put 

(6.30)  Q°2  =  (D  ^^  D  ^^  q) 

P    P     P=l 
so  that 
(6.30a)   qll   =  X[p°2]  5.^  +  tx(p°2  ^  p°2)  ^.  qo2  ^ 

By  (6.28) 

/  o2   „o2v   ^  ,  o2   „o2s    _  ^ 

and  hence,    since  q°      -  Q°      is   symmetric, 

(6.31)  (Q°^  -  Q^^)^!  +   (Q^3   -  Q32)^2  =  ^ 
By  (6.30) 

^^ik  -  ^ki^^ik  =   t<^°^(^''-  ^)1 

=   (°  ol°  01^^   U(P)-D  ol^  oA--U(P)Jp=l 
P       P  P       P  ^ 


=   U3(p°l,p°\r    -  <^-) 


rrt: 


P      .r.     '■-■■f     V- 


,*/»  k  fj.  ••,  /     '•- 


■•'i-  V 
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Since  E.    -  ?;"  is  skew  symmetric,  the  terms  with  A,B,C,  in  the 
expression  O.l^b)  for  U^{p°  ,p°  ,f:  -  E'^)   make  no  contribution, 
and  we  find 

(^?k  -  ^kf^^ik  =  2X[p°l][p°l(^^  E)]   +2ix[p°V^(c''-  a) 
identically  in  f .   In  particular 

3i  '  13  ~   ^   ^^^31  "  ^13^   Sixlp^.  p.^  -  pj^.  Pj^i 
For  i=l,2  we  obtain  then,  using  (6.21;) ,  (6.13) 
Q3I  -  Q?^  =  ^X[p°l]v^,  .2,(v^,(pJ^.p°l) 


=  2(q°|-v^^  +  q°lv^2  "^  ^3^,1^ 


Hence  by  ( 6.15) , (6.17) 


Q' 


o2   ^o2 


31  -  Q13  =  2(4^22  ^,1  -  ^12  ^,2^ 


Q°|  .og  =  2(-4,2i  \1^<^,11\2)- 

Substituting  these  expressions  into  (6. 31)  we  obtain  the  do- 
sired  differential  equation  for  v: 

(6.32)    v^-L^  4:^22  -  2v^^2  ^11  ^^'22  ^'11  =  0- 

Unless  4  vanishes  identically  the  differential  equations 

(6.23) ,(6.32)  describe  the  behavior  of  tangential  stresses  and 

2 
normal  displacements  to  lowest  order.   Up  to  terms  of  order  h 

the  stress  matrix  tis  given  by  hq   and  thus  by  ( 6.15) ,(6.17) 
^11^-^^,22'  '-12  =  ^21  ''-  ^il2»  '^22^'^^,11 

H3  =  ^31  ^"  ° 

where  4  Is  a  solution  of  the  bl-harmonic  equation.   Thus  the 

2 
stresses  normal  to  the  surface  are  of  order  h  while  the  tan- 
gential stresses  are  of  order  h,  unless  the  second  derivatives 
of  4  vanish.   The  differential  equation  for  the  tangential 
stresses  does  not  involve  the  vertical  deflection, since  here 
for  i=l,2 


S 


'f.      s 


ofi-r 


'   h 


r.- 
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t 

.   -^  =  P31  =  hP3i  +  O(h^)  =  hv^^  +  O(h^) 
the  vertical  deflection  is  given  by 
(6.33)    x^  =  hv{x^,x^)   +  O(h^)  , 

where  v  satisfies  equation  (6.32),  Thus  up  to  terms  of  higher 
order    000 

•^v'  >*^-J-'  ..^v' 

We  next  consider  the  case  where  all  stress  components  q.. 
are  of  the  order  h  i.e.  where 

(6.31;)    q°-^  =  0. 

In  this  case  the  differential  equation  (6.32)  collapses  and 
higher  order  terms  In  the  expansion  have  to  be  considered. 
It  follows  from  (6.I8)  that 


r^ol  _  ^ol  _  „ol  _  ^ol  .  ^ol  _  ^ 
Pll  -  P22  -  P33  ■  P12  ""  P21  -  °- 


Then  also 


Hence 


^ol   _   ^ol   _  „ol   _  ^ 
Pl2,l  -  -  P21,l  -  Pll, 2  -  ^' 

Pl2,2  "  P21,2  "  "  P22,l  "  ^' 


p,  P  =  -Pp-i  -   const. 


By  a  suitable  rigid  motion  in  the  x-jXo  -plane  one  can  bring 
about  that  this  constant  vanishes.  That  leaves  only  the  com- 
ponents 

_ol  „ol  ^ol  _ol 
31*   32'   13*  23 

in  p°  which  are  described  by  the  single  fimction  v(x,,Xp) 
through  equations  (6.2[j.),  where  by  (6.33)  hv(x,,X2)  is  essen- 
tially the  normal  displacement  of  the  plate.  We  notice  that 
now  p°   is  skew  symmetric: 

(6.3^)    P°^  =  -  P^^''  . 


t-F,A 


it.    : 
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Then  by  (6.3O) , (3 .l5b) 

Qll=    [0°2^"]  =  (D,iD^,  D  U(p)) 

P    P    ^      P=l 

=  u3(p<^\p°l,,) 

(6.36)  =  x[r][p°iV^]  +^I  [(^  +  ^'')  P°^V^] 

Substituting  the  specific  values  of  the  p?.  in  terms  of  v,  we 
find 

^?k  =  2Mv^,2^  ^/)    5^^  -.  2^   v^,  v^^  for  i=k=l,2 

(6.37)  Q°3  =  Q3f  =  0  for  1=1,2 

Q°^  =  (2X  +  2tx)(v^^2+  V  2^). 

Equations  (6.37)  show  that  ^°   is  symmetric.   It  follows 

from  (6.30a)  that  also  the  matrix  q°  is  symmetric,   ^-'e  have 

then  from  (6.28) 

^o2    .  „o2   _  p,   _o2    .  ^o2   _  ^   ^o2  _  „o2 
^11,1  ""  ^12,2  -  ^'    ^21,1  ^  "^22,2  -  °'   ^12  "  ^21 

It  follows  from  these  equations  that  there  exists  an  Airy 
function  i)f  =  ^{/(x^jXo)  such  that 

i^   oA\    ^o2  _  .    _   ^o2  _  „o2  _  ,     „o2  _  , 

(6.38)  q^^  ~    ^^22  -  '  ^12  -  ^121  -  1',12»  *^22  "  ^ll' 

Prom  (6.28)  and  the  symmetry  of  q°  we  have  for  the  remaining 
components  of  q 

(6.39)  q°3  =  q3^  =  0   for  1=1,2,3. 

We  obtain  a  differential  equation  for  \^  by  applying  iden- 
tity (6.22)  to 

E   =  p°^  ?=  (D  ^2  ^^P))   .  =  q°^  -  Q°2. 

p        p=l 

The  right-hand  side  in  (6.22)  vanishes  because  of  the  compat- 
ibility conditions 

/  ^  I  /^\     o2   _  ^o2 
(6.1;0)    p^2,l  -  Pil,2 

for  p°^;  (see  (6.5)).   Substituting  for  the  q^^  and  Q°^  their 
val\aes  from  (  6.38) ,(  6.39) ,(  6.37)  we  obtain 


•i,V 


t .  . , 


,     f 


<^■ 


63 
A  second  equation  between  \\i    and  v  is  derived  from  the  fact 
that  the  relations  connecting  q°-^  and  p°  have  to  be  compat- 
ible. We  find  from  (6. 6), (6. 7)  for  y=3 


1    '"11      '  "li,i     ^leifd.  "11,1   l<i,<i 

for  1=1,2,3. 


A  combination  of  these  equations  is 

'^31  -  ^W'.l  *   '03I  -  ^l3',2 

This  relation  will  yield  the  desired  differential  equation. 
We  find  as  in  ( 5.1l4-a,b,c,d,e,f ) 

(6.1,3a)    q°^  =  (D  01  °  01  °  01  ^  -^  3D  ^-L  D  ^2  q  +  D    q) 

P    P    P         P    P        P     P=l 

(6.1,3b)    q21  =  (D  1^  D  ^^  D  -L^  q 

P    P    P 

■^  °  01  °  20  ^  ^  2D  ^^  D  13^  q+D  21  q) 
p    p         p    p      p     p=l 

or  for  any  matrix  ? 

(6.l,l,a)    [q°-^  e"]    =  U^(p°\p°\p°\r) 
+  3u3(p°\p°^£)  +  U^p^^r) 

(6.l,l,b)    [q2V-^=-]  =  U^(pl^pl^p°^n  +  U3(p°\p2°,?) 

+2u3(p^°,p^^,r)  +  u2(p21,r). 

Since  only  combinations  q?^  -  q?^  enter  (6.I4.2)  we  only  need  to 
determine 

[(q°3  -  q°3")r'^]  =  [q°3(^'^-  ="  )  3 

=  U^(p°^p^\p°\g  -f^=-)  .  3u3(p°l,p°2,r  -  f*) 
+  U2(p°3,^  .  ^-^^-). 

Using  the  expressions  (3.l5a,b,c)  for  ir,U-^,ir"  we  see  that  only 
terms  with  coefficients  A  or  (j,  make  contributions,  since  both 
p°   and  ^-f'"'  are  skew  symmetric.   Ve  obtain 


,  r:-. 


»•.  ;  » 


;^     f-?-. 


•...L.1 


^:  o 
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Here  p°   is  the  matrix  given  by 

/O    0    -V  A 

(6.i^5b)   p°^  =  '  0    0    -V  ,1 

I  '  / 

Taking  for  ^   the  special  matrix 

/   0  0  0  \ 

C  =  I   0       0       0  ! 


1 


\    1  0  0  / 

relation   (6.14.5)    reduces   to 
(6.l;6)  q°l   -  q^^  =   +  12(  A+ti)(  v^j^2+v^/)v^  ^^  +  6{\+[i)[v°^]v  ^^ 

The  p?.    are    connect.^d  with  the  known  q.,    by   relations    (6.30a) 

which  yield 

o2        r  -,\    ,    rt    \      o2    .    ^o2 
q        =    (3\  +  2p,)    p        +   Q 

/,o2  _   ,    _o2    .    o  ^o2   _^  ,^o2 
^22        '^  P        "^  2[iP22   +  ^22 

^o2  _      /    o2    .      o2x    .    p,o2 
qi2  =   ^^(Pi2   -^  P21)    ^   ^12 

Hence,    using  the   values   of  the   q?^^,    Q°j^  given  by   (6.37),    (6.38), 

(6.39) 

[p°2]   =   .  2(v  ^2+  V     2)   + 


,1    •    \2    '        3A+211"'' 
2m.  P^I  =   ^l/^ii   -  IXT^    ^*   -  2ix  v^2^ 

ix(p°|  +  P^i)    =   *^i2   -  2t^^,iV,2 
Substituting   these   values    into   (6.I46)    we   obtain 

^31   -   ^13  =  ^^^^,22   ^,1   -^12  ^,2^- 


vl" 


..  V  •. 
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Similarly,  it  follows  that 

^32  -  ^23  =  ^^^11  \2-^l2  \l)- 
This  enables  us  finally  to  express  the  left  hand  side  of  (6.I4.2) 
in  terms  of  \j;  and  v: 

(6.^7)       (q3l-q?^),iMq°|-q°3)^2 

"  3<*,22  ",11  -  2*,i2  v_j2  +  >|r_^j^  v^^g' 

It  remains  to  carry  out  a  similar  reduction  for  the  right 
hand  side  of  (6.14.2).   The  matrix  q   is  determined  by  (6.l4l4.b). 
Since  by  (6.26 ),( 6.27)  p^°   and  p^°  vanish  the  terms  with 
U^  and  U-^  in  the  expression  (6.Li.i;b)  make  no  contribution  and 
we  are  left  with 


or 


qfj=  Up21je,^  +  ^(pg.p2];). 


21 

Thus  q   is  symmetric.   Moreover  by  (6,[|.)  for  a=p=l 

^21  _   „11     „11 
^i3  "  "  ^il,l  '  ^i2,2' 

Hence 

(6.1.8)  34]    -    q^l  =  2qfl=  -2(q]:l,i.  qj^,^). 

By  (5.1i;e)   for  x?-^   =  0 

(6.1.9)  q^^  =  (D  ^3^  q)     =  X[p^-]  +  ix(p^^  +  p^^"^ 

P     P=l 
Here  from  (6.1,)  for  0=0,  p=l 

(6. 50)  qj^  =  q°];^^  -  q°^^2  =  ^    ^°^^  i=l,2,3, 
since  q°  =  0.   In  particular  for  i=3 

(6. 51)  0  =  q^^  =  X[p^^]  +  2,tp^^  . 
On  the  other  hand  by  (  6.5)  5  (  6.i|5a) 

(6.52)  PiJ  =  P?3,i,  =  -\ik   for  i,k=l,2 

[p^^]  =-Av  +  p^^ 


^^4 


"0 


^-,  J. 
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and  thus  by  (6.51) 

(6.53)   [p^^]=  xJI^Av. 

Combinlns  (6.14.9)  with  (6.52) ,( 6.53)  we  get 

^Ik  =  -lis  ^k^^  -2|xv^,j^  for  i,k=l,2. 
Substituting  these  values  for  the  qfr;  Into  (6.14.8)  we  find  that 
-,„21   ^21   +8u,(X+a)       p   «_i  o 

and  hence  that 

(6.5I1)     (3qi3  -  q3i)^i  +(3q23  -  o,^^^  ^2.   '   -^-^^^   ^- 
Combining  (6.514-) ,( 6.I4.7) ,  (6.I4.2)  we  obtain  the  desired  second 
differential  equation  connecting  v  and  '^'. 

(6.55)  ^^   V  =  |i^Jfi}  (*^22  v_ii  -  2*^12  v_,2  .♦_,,  v^22) 

Up  to  terms  of  higher  order 
w  =  hv(x,  ,Xp) 

represents  the  normal  displacement  of  the  plate,  while 
P        ?  2  2 

'/       _  "Ci       oZ  _   h"^,     -,<.   _   h^,     ;^  ~^\, 

^11  ~  "2  ^11  ~  2  ^,22'  42  ~  '  "2  '♦',12'  '  22  ~  :?  ^,11 
represent  the  tangential  components  of  the  stress  in  the  middle 
surface.   These  quantities  represent  at  the  same  time  the  aver- 
age tangential  stresses  throughout  the  thickness  of  the  plate, 
since  q  °  =  0  and  the  average  of  q  x^h  vanishes.   Equation  (6.5) 
can  be  written  In  the  form 

(6.56)  A  w  =  -^^t^^   w^ii  +  2ti2  ^,12  ^  '^22   ^,22^ 
where 


»'  =  t^>'' 


is  thfe  plete  stiffness  constant,  ^'e  Introduce  the  Airy  stress 
function  2 

so  that 


"^11  "^,22  '   ^12  "^   ^21  "  "-^,12  '  '^22  "  ^,11  • 


-  \ 


iC-..   '■<■ 


07    0"J 
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Then  equation  (6,56)  takes  the  form 

By  (6.1^.1)  w3  have  as  second  equation  between  X  and  w 

(6.58)    <^  ^  =  £(>!  -,2  w  2-1  -w  22  ^  11^ 
where 

E  =  M-(3^+2|i) 

is  Young's  modulus,   Equations  ( 6.57) , (6.58)  represent  the  fam- 
iliar V.  Karman-Poppl  system  for  buckling  of  plates. 


ai- 


References 

[1]  V.V.  Novozhilov:   Foiindatlons  of  the  nonlinear  theory  of 
elasticity.  (Graylock  press,  1939). 

[2]   P.D.  Murnaghan:   Finite  deformation  of  an  elastic  solid, 
(viiey  Sons,  1951). 

[3]  W.Z.  Chlen:  The  intrinsic  theory  of  thin  shells  and 
plates.  Quart.  Appl.  Math.  l(19i4-3)  pp. 297-327, 
II,  III  ibid.  V.2,  19kh,   PP.U3-59,  120-135. 

[I).]   J.L.  Synge  and  W.Z.  Chien:   The  intrinsic  theory  of  elastic 
shells  and  plates.   Applied  Mathematics,  Th.v. 
Karman  anniversary  volume  (I9I4-I)  pp.  103-120. 

[53  E.  Reissner:   On  a  variational  theorem  for  finite  elastic 
deformations.   J.  Math.  Physics  32(1953)  pp. 129-135. 

[6]  R.S.  Rlvlin:  Lectures  on  solid  mechanics.  Seminar  in 
Applied  Mathematics,  Boulder,  Colorado,  1957  (in 
preparation). 


»•,  -K.  . 


r      .'. 


.i^t 


£J 


»  >^<.-*  "^  - 


.1:    r 


I-  •  r.  - 


•■-       ryr,^ 


J  .  ?f  »  .J 


fT    .?      fM 


DISTRIBUTION 
CONTRACT   NO. 

Commander 

Armed  Services   Technical 

Information  Agency 
ATTN:      TIPDR 
Arlington  Hall   Station 
Arlington  12,    Virginia  10 

CoramandGr 

Air  Force  Ballistic  Missile 

Division 
Air  Research  +  Development 

Command 
P.    0.    Box  262 
Inglewood,    California  1 

Commander 

Air  Force  Office  of  Scientific 

Research 
Air  Research  +  Development 

Command 
ATTN:  Mechanics  Division 
Washington  25,  D.  C.  2 

Director  of  Advanced  Studies 
Air  Force  Office  of  Scientific 
Research 

Box  2035 

Pasadena  2,    California  1 

Commander 

Air  Force  Armament  Center 

ATTN:  ACOP 

Eglin  Air  Force  Base,  Florida  1 

Commander 

Arnold  Engineering  Development 

Center 
ATTN:   Library 
P.  0.  Box  162 
Tullahoma,  Tennessee  2 

Commander 

Air  Force  Cambridge  Research 

Center 
L.   G.    Hanscom  Field 
ATTN:      CRQ3T-2 
Bedford,   Mass.  1 

Office   of  Technical   Services 
Department  of  Commerce 
Washington  25.    D.    C.  1 


OF   REPORTS 

AF~WTT3HT-i6i 

Commander 

Air  Force  Fllgth  Center 
Edwards  Air  Force  Base 
California  1 

Commander 

Air  Force  Special  Weapons  Center 
■  Kirtland  Air  Force  Base 
New  Mexico  1 

Commander 

Air  Research  +  Development  Command 

ATTN:   Director  of  Research 

P.  0.  Box  1395 

Baltimore,   Maryland  1 

Commander 

Hq.    Air  Force  Missiles   Development 

Center 
Holloman  Air  Force  Base 
New  Mexico  1 


Commander 

Air  Force 

Ml 

.ssile 

Test 

Center 

Patrick  A: 

Lr 

Force 

Base 

Florida 

Commander  5 

Wright  Development  Center 
ATTff:      Aircraft  Laboratory    (1) 
Aeronautical  Research 

Laboratory   (1) 
Director  of  Weapons 

Systems  Operations    (1) 
WCOSI-3    (1) 

Materials  Laboratory    (1) 
Wright-Patterson  Air  Force  Base, 
Ohio 

Commanding  Officer 

Maxwell   Air  Force  Base 

ATTN:   Air  University  Library 

Maxwell  Air  Force  Base,    Alabama 


A  c.  r.mo  Ci 


"-JO  rui'.';  0 


X'jB 


J    .  1 


*• 


Conunander 

Air  Materiel   Cornmand 
Wright-Patterson  Air  Force  Base 
Ohio  ] 


Office  of  the 
Department   of 
ATTN:      ORDTB 
Washington  2S, 


Chief  of 
the  Army 

D.    C. 


Ordnance 


Director 

Ballistics   Research  Lab. 
Aberdeen  Proving  Ground 
Aberdeen,   Maryland  1 

Commanding  Officer 

Office  of  Ordnance  Research 

Duke   Station 

Durham,    North  Carolina  1 

Commanding  Officer 

Redstone   Arsenal 

Huntsville,    Alabama  1 

Commanding  Officer 

White  Sands  Proving  Ground 

Las   Cruces,    New  Mexico  1 

Commander 

Air  Force  Office  of  Scientific 

Research 
Air  Research  +  Development 

Command 
ATTN:   Directorate  of  Research 

Information 
Washington  2S,    D.  C.  2 


Chief  of  Naval  Research 
Department  of  the  Navy 
ATTN:  Mechanics  Branch 

Air  Branch  (1) 
Washington  2S,    D.  C. 


(1) 


Director 

U.S.  Naval  Research  Laboratory 
ATTN:  Mechanics  Division 
Washington  25,  D.  C. 


Chief,  Bureau  of  Aeronautics 
Department  of  the  Navy 
ATTN:   Research  Division 
Washington  ZS,    D.  C.  1 

Chief, 

Bureau  of  Ordnance 

Department  of  the  Navy 

ATTN:   Research  +  Development  Div. 

Washington  25>  D.  C.  1 

Commander 

Naval  Ordnance  Laboratory 

White  Oak 

Silver  Spring,  Maryland        1 

Commander 

Naval  Ordnance  Test  Station 

Inyo kern 

China  Lake,  California         1 

Superintendent 

Naval  Postgraduate  School 

Monterey,  California  1 

Commander 

David  Taylor  Model  Basin 

Washington  7,  D.  C.  1 

Director 

Aeronautical  Research 

National  Advisory  Committee  for 

Aeronautics 
1^12  H  Street,  N.W. 
Washington  25,  D.  C.  1 

Director 

High  Speed  Plight  Research 

Station 
National  Advisory  Committee  for 

Seronautics 
Edwards  Air  Force  Base, 
California  l 

Commander 

European  Office 

Air  Research  +  Development  Command 

k-1   rue  Cantersteen 

Shell  Building 

Brussels,  Belgium  l 


.c.-f... 


-3- 


Dlrector 

Ames  Aeronautical  Laboratory 

National  Advisory  Committee 

for  Aeronautics 
Langley  Field,  Virginia        1 

Director 

Lewis   Flight  Propulsion  Laboratory 

National  Advisory  Committee  for 

Aeronautics 
Cleveland  Airport 
Cleveland,  Ohio  1 

Director 

Langley  Aeronautical  Laboratory 

National  Advisory  Committee 

fpr  Aeronautics 
Langley  Field,  Virginia        1 

Director 

National  Bureau  of  Standards 

Washington  25,  D.  C.  1 

U.S.  Atomic  Energy  Commission 
Technical  Information  Service 
1901  Constitution  Ave. 
Washington  25,  D.  C.  1 

Office  of  Technical  Services 
Department  of  Commerce 
Washington  25,  D.  C.  1 

Applied  Mechanics  Reviews 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio  6,  Texas  2 


Aeronautical  Engineering  Review 

2  East  61+th  Street 

New  York  21,  New  York  2 

Professor  George  Herrmann 
Dept.  of  Civil  Engineering  and 

Engineering  Mechanics 
Columbia  University 
New  York  27,  N.  Y. 


Brown  University 

Div.  of  Engineering 

Providence  12,  Rhode  Island     1 

University  of  California 
Low  Pressures  Research  Project 
Inst,  of  Engineering  Research 
Berkeley  l{.,    California         1 

California  Inst,  of  Technology 

Jet  Propulsion  Laboratory 

I|.800  Oak  Grove  Drive 

Pasadena  3,  California         1 

California  Inst,  of  Technology 
Guggenheim  Aeronautical  Lab. 
Pasadena  k,    California         1 

Cornell  University 

Graduate  School  of  Aeronautical 

Engineering 
Ithaca,  New  York  l 

John  Crerar  Library 

Chicago   1,    Illinois  1 

Harvard  University 

Dept.  of  Engineering  Sciences 

Cambridge  38,  Mass.  1 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
8621  Georgia  Ave. 
Silver  Spring,  Maryland        1 

Dr.  Francis  H.  Clauser 
Dept.  of  Aeronautics 
Johns  Hopkins  University 
Baltimore  18,  Maryland         1 

Marquardt  Aircraft  Corporation 

7801  Havenhurst 

Van  Nuys,  California  1 

AVCO  Research  Laboratory 
ATTN:   Chief,  Technical  Library 
2385  Revere  Beach  Parkway 
Everett  i^9,  Mass.  1 


-k- 


AVCO  Manufacturing  Company 

155  Sniffins  Lane 

Stratford,  Conn.  1 

Bell  Aircraft  Corporation 

ATTN:   J.S.  Isenberg 

Box  1 

Buffalo  5,  New  York  1 

Boeing  Airplane  Company 

Box  3107 

Seattle  II4.,  Washington        1 

Chance-Vought   Aircraft  Corp. 
Dallas,    Texas  1 


Convair 

Port  Worth  Division 

Port  Worth,  Texas 

Convair 

Pomona,  California 


1 
1 


Convair 

San  Diego  Division 

ATTN:   Chief,  Applied  Research 

San  Diego,  California         1 

Cornell  Aeronautical  Lab.,  Inc. 

U-k-SS   Genesee  Street 

Buffalo  21,  New  York  1 

Douglas  Aircraft  Company,  Inc. 

3000  Ocean  Park  Boulevard 

Santa  Monica,  California       1 

Pairchild  Engine  and  Aircraft 

Company 
Guided  Missiles  Division 
Wyandanch 
Long  Island,  New  York  1 

General  Electric  Company 
Aircraft  Gas  Turbine  Division 
Cincinnati  15,  Ohio  1 

General  Electric  Company 
Special  Defense  Products  Div. 
3198  Chestnut  Street 
Philadelphia  k,    Pa.  1 


The  Glenn  L.  Martin  Company 
Baltimore  3,  Maryland         1 

Armour  Research  Poimdation 
Illinois  Institute  of 

Technology 
Chicago,  Illinois  1 

Qvimman   Aircraft  Engineering  Corp. 

Bethpage 

Long  Island,  New  York         1 

Dr.  A.  E.  Puckett 

Research  +  Development  Labs. 

Hughes  Aircraft  Company 

Culver  City,  California        1 

Lockheed  Aircraft  Corporation 

Box  551 

Burbank,  California  1 

Lockheed  Aircraft  Missile 

Systems  Division 
Van  Nuys,  California  1 

McDonnell  Aircraft  Corporation 

Box  516 

St.  Louis  3,  Missouri  1 

University  of  Maryland 
Institute  for  Fluid  Dynamics 

and  Applied  Mathematics 
College  Park,  Maryland         1 

Naval  Supersonic  Laboratory 
Mass.  Institute  of  Technology 
Cambridge  39,  Mass.  1 

Library    (33-30L).) 

Mass.    Inst,    of  Technology 

Cambridge   39,    Mass.  x 

Library 

Institute  of  Aeronautical  Sciences 

2  East  6Uth  Street 

New  York  21,  New  York  i 

University  of  Michigan 
Department  of  Aeronautical 

Engineering 
East  Engineering  Building 
Ann  Arbor,  Michigan  1 


;:r 


ri"^: 


.<-. . 


>..»... 


-5- 


University  of  Minnesota 
Engineering  Library 
Institute  of  Technology 
Minneapolis,  Minnesota        1 

University  of  Minnesota 
Rosemount  Research  Center 
Rosemount,  Minnesota  1 

New  York  University 

Dept.  of  Aeronautical  Engineering 

University  Heights 

New  York  53,  New  York  1 

Ohio  State  University  Research 

Foundation 
Columbus  10,  Ohio  1 

Polytechnic  Institute  of 

Brooklyn 
Dept.  of  Aeronautical  Engineering 

and  Applied  Mechanics 
99  Livingston  Street 
Brooklyn  1,  New  York  1 

Polytechnic  Institute  of 

Brooklyn 
Aerodynamics  Library 
527  Atlantic  Avenue 
Freeport,  Nex-j  York  1 

Professor  S.  Bogdonoff 

The  James  Forrestal  Research 

Center 
Princeton  University 
Princeton,  New  Jersey         1 

Rensselaer  Polytechnic  Inst. 
Dept.  of  Aeronautical  Engineering 
Troy,  New  York  1 

University  of  Southern  California 

Engineering  Center 

3518  University  Avenue 

Los  Angeles  7,  California      1 

Stanford  University 
Guggenheim  Aeronautical  Lab. 
Stanford,  California  1 


The  University  of  Texas 

Defense  Research  Laboratory 

P.  0.  Box  8029 

Austin  12,  Texas  1 

Linda  Hall  Library 
ATTN:   Mr  T.  Gillies 

Document  Division 
5109  Cherry  Street 
Kansas  City  10,  Missouri       1 

Dr.  H.  G.  Lew 

Dept.  of  Aerona\itical  Engineering 
The  Pennsylvania  State  University 
University  Park,  Pa.  1 

Aerojet  Engineering  Corporation 

6352  N.  Irwlndale  Avenue 

Box  296 

Azusa,  California  1 

North  American  Aviation,  Inc. 
Aerophyslcs  Department 
1221/4-  Lakewood  Boulevard 
Downey,  California  1 

Northrop  Aircraft,  Inc. 
Hawthorne,  California         1 

The  Ramo-Woolrldge  Corporation 

ATTN:   Chief  Librarian 

5730  Arbor  Vltae 

Los  Angeles  k5,    California     1 

Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California       1 

Reed  Research,  Inc. 

10[).8  Potomac  Street,  N.W. 

Washington  7,  D.  C.  1 

Republic  Aviation  Corporation 

Farming  dale 

Long  Island,  New  York         l 

United  Aircraft  Corporation 

Research  Department 

362  Main  Street 

East  Hartford  8,  Conn.         1 


) 


iilo't.'rljBl'i   :; 


•6- 


Westinghouse  Electric  Corporation 
Aviation  Gas  Turbine  Division 
Lester  Branch  Post  Office 
Philadelphia  13,  Pa.  1 

University  of  Washington 

Dept.  of  Aeronautical  Engineering 

Seattle,  Washington  1 

Stanford  University 

Dept.  of  Aeronautical  Engineering 

Stanford,  California  1 

University  of  Toronto 
Institute  of  Aerophysics 
Toronto  5,  Canada  1 

Dr.  Leon  Trilling 
Mass.  Inst,  of  Technology 
Fluid  Dynamics  Research  Group 
Cambridge  39,  Mass. 


1 


i 


) 


y 


..70 


CA 

NYU 

IMM- 

250      John. 

On  finite  deformations  of  an 
elastic  isotropic  material. 


i^ 


NYU 
2p0 


c.  1 


John 


AUTHOR 

On  finite  deformations  of 

TITLE 

an  elastic  isotropic  materia] 


DATE  DU  E 


^ 


BORROWER  S  N^E 


/4^C>(c' 


^ov  .s-*^ 


N.  Y.  U.  Institute  of 
Mathematical  Sciences 

25  Waverly  Place 
New  York  3,  N.  Y. 


